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SUMMARY 


This paper presents a method whereby the forced vibration 
characteristics of a system formed by interconnecting a group 
of subsystems can be predicted in terms of the interconnection 
constants and properties of such subsystems which are readily 
and accurately determinable by experiment. These properties 
are the generalized dynamic admittances—motions of system 
components due to application of unit oscillating forces. The 
method was devised to enable selection of a means of connecting 
two structures together so that excessive vibration due to im- 
posed forces can be prevented, under the conditions that: (1) the 
structures already exist and cannot be modified; (2) the inter- 
connection constants are the only quantities over which selection 
can be exercised; (3) the structures are, in general, so complicated 
that formulation of their equivalent dynamical systems and 
analysis thereof by conventional methods is impracticable or 
impossible. 

Elementary tensor and matrix algebra comprise the mathe- 
matical tools. Application of the theory to aircraft power-plant 
installation problems is discussed, but extensions to other prob- 
lems and fields are possible. 


INTRODUCTION 


HE BEHAVIOR of a vibrating continuous system 
BD can be specified with sufficient accuracy by the 
equations applicable to a multiple degree-of-freedom 
lumped parameter system. Attention is restricted 
herein to systems excited by sinusoidal forces of 
identical frequency for which the differential equations 
admit the solution: 


[(Q) = [oullad (1) 


Here Q and g are generalized complex forces and dis- 
placements, respectively. The term “impedance matrix’’ 
will be applied to [p.;) whose typical element is given by 


Pig = Cig — Aygo? + Oi (2) 
Solving for the displacements: 
[ad = [oul 1Q] = [Yul [21 (3) 
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[Yi] ts called the admittance matrix. Values of its 
elements, which are functions of frequency, can be found 
experimentally by applying a unit oscillating force Q; 
to each coordinate q; and measuring the resulting ampli- 
tudes and phases of the coordinates q,. Thus the ratios 
q:/Q; are the admittances, classified as direct or cross 
according asi = jandi# ij. 6;;(= argYj,) is the phase 
angle by which q; leads Q;._ Resonances of the system 
are indicated by maxima of the Y;j,’s. Further 
properties of admittances per se are well known and 
will not be discussed here.': * 

In practical studies of complicated systems requiring 
a large number of coordinates for their specification, 
interest is usually limited to a reasonable fraction of . 
the total coordinates. For such cases, admittance 
methods become exceptionally attractive, permitting 
an arbitrary number of displacements to be found in 
terms of the known applied group of forces when the 
appropriate admittances have been experimentally 
determined. 

The purpose of this paper is to reveal the charac- 
teristics of the system formed by combining a group of 
subsystems whose admittance properties are known. 
This is the problem imposed, for example, when a 
machine, the operation of which produces oscillating 
forces, is to be attached to a complex structure. Criteria 
relating to the mounting of engines in aircraft are par- 
ticularly sought, and obviously no loss of gener- 
ality will ensue if illustrations from this field are 
used. 


CLASSIFICATION OF DYNAMICAL SYSTEMS 


Since the development and application of the follow- 
ing theory depend to a large extent upon recognition of 
the type of system under consideration, some space 
will be devoted to a discussion of various characteristic 
systems and the explicit qualitative form of their 
impedance matrices. 

Eqs. (1), (2), and (3) are to be interpreted hereafter 
as follows: the displacements or configuration of a set 
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of n rigid bodies will be given by simultaneous values 
of g;, where each gq, is a compound coordinate, or a set 
of six absolute displacements of the body 7 referred 
to a rectangular coordinate system, the origin of which 
is located at the equilibrium position of the center of 
gravity of body 7. The m coordinate systems are 
established so that their corresponding axes are paral- 
lel and have positive extensions in the same di- 
rections. A typical system is illustrated in Fig. 1, 
for which q, = {x1, Wy, 2a, i Wy, 6} and gq = 
{ x2, Ya, Ze, G2, V2, Oo}. 
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Dynamical systems requiring more than one co- 
ordinate for the specification of each body will be 
called ‘compound.’ The quantities g, p, and Q are 
then compound tensors. Systems in which the bodies 
are constrained to move in one and the same direction 
only or.in which the interconnections are such as to 
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decouple some one mode from the others are to be 
called “simple,” and g, p, and Q are then ordinary 
tensors with scalar elements. 

Both simple and compound systems are subject to 
one of two classifications, depending on the manner 
in which the various bodies are joined. If each body 
except the first and last is elastically connected with 
only two others, the combination is called a ‘‘chain 
system.” If there are more interconnections, no 
special identification will be assigned. The form of the 
impedance tensor reveals at a glance the classifications 
described. 

Table 1 contains an illustration of each classification 
and a qualitative or schematic representation of the 
appropriate impedance tensor. Blank elements denote 
zero va'ues, and, for the coordinate system chosen, 
only diagonal elements of the diagonal tensors have 
the form py; = C4 — a@;w? + jbiw. The other nonzero 
elements are pi; = ¢;; + jbi. Dashpots contributing 
b terms are not shown on the diagrams. Absence of 
certain elements of the diagonal tensors in the com- 
pound systems illustrated indicates decoupling of the 
corresponding coordinates due to symmetrical spring 
configurations. It is helpful to recall that in each case 
the dynamic tensor is the sum of the stiffness tensor c, 
jw times the damping tensor, and —w? times the inertia 
tensor. 

For the chosen coordinate system there is only 
static coupling, and the elements of the stiffness 
tensor can be found from 0?V/0g,0q;, where V is the 
potential energy. Similarly, the dissipation function 
F enables the elements of the damping tensor to be 
calculated by 0?F/ 04,04). 

Investigation of the dynamic characteristics of the 
configuration formed by connecting together two sub- 
systems (in terms of the known subsystem properties 
and the manner of interconnection) can now begin. 
It will be found that a development applicable to com- 
pound systems, involving the use of tensor principles, 
introduces mathematical advantages of generality and 
unity which overbalance the physical complexities. 
In fact, certain properties of simple systems which 
cannot be revealed by matrix methods are easily found 
as a special case from the general treatment and will 
be given later. 


ADMITTANCES OF COMBINED COMPOUND SYSTEMS 


In the following presentation, compound tensors 
are expressed mainly in matrix notation, identifying 
subscripts being employed to prevent ambiguous 
assignment of indices. Frequent use is made of direct 
notation, but in no case does a subscript denote co- 
variant properties—that is, true index notation is not 
used. 

‘Consider two systems, S; and S:, with compound im- 
pedance tensors: 
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 pirpie . + +Pin 
P21P22- - - P2n 
an = 
| PniPn2- - -Pnn 
, (4) 
P(nt+1)(n+D P(n+1)(n+2)- - -P(nt+1) (n+p) 
P(n+2)(n+1) P(n+2)(n+2)- - - P(n+2)(n+p) 
a = 
LP (n+p) (nt+1) P(n+p)(nt+2)- - -P(ntp) (n+p) 


A typical element of p; or p: hereafter will be under- 
stood to have the form: 


Ixlylz|¢lylel 
x [Aa] Ay 
Y |x |Ay| * e e 
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The complex elements thereof are given by Eq. (2). 
If the systems are connected in an arbitrary way, the 
combined impedance tensor, p*, will have elements 
that are not, in general, simple functions of the original 
subsystem parameters. However, let them be con- 
nected by means of a known arrangement of springs 
and dashpots between the mth body of S; and the 
(n + 1) body of S:. Introducing compound coordi- 
nates g, and gg, the connecting spring-dashpot im- 
pedance tensor, pz, can be specified: 


[FT Ge | 
fra* lel ALL 
%\ Le KX’ (5) 
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The collection of systems S, and S; and the connecting 
elements can be represented by the tensor, p, before the 
interconnection is made. 
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Interconnection is now made by merging the coordinate 
fq With g, and gg with g,+:1 and is defined by the 
equations 


q = Tq* 


where g* are the new coordinates, and T is the singular 
transformation tensor. 
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The impedance of the combined system is found from 


p* = TT 
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In expanded form: 
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In practice, the systems S, and S; are too complicated 
to permit determination of reasonable values of the 
impedance elements. The representation of compli- 
cated continuous structures, such as aircraft fuselages 
and wings, by a sufficiently large number of concen- 
trated masses, in order to study the behavior of various 
small components, is often an unsatisfactory com- 
promise between accuracy and effort in which the 
parameters cannot but be uncertainly selected. How- 
ever, any element of the admittance tensor, the re- 
ciprocal of the system impedance, can be found experi- 
mentally with great accuracy and little effort. 

The mathematical problem, then, is to find the 
combined system admittance tensor, Y*, in terms of 
the known tensors Y,, Yo, K, K’, and L. The fact 
that all the elements of Y; and Y2 need not be known 
creates no difficulty. 

If a square matrix, a, and its reciprocal, 8, be parti- 
tioned according to the scheme: 


aa | ee] ond ect eg = | Pubs 
@21! 29 Bor! Boe 


it can be shown*® that the elements of 8 are given by 
the alternative forms 








Bu = ana! + an! ang"(a22 — aera !ay2) ~ eta! = (ai, — ag*a22~!*a1) ~? 

Bis = — a ar2(a%22 — aera an2)~! = — (arr — @12"Q22~! a1) ~ er" O22" (9) 
Bor = —(a2 — aera ay) agra! = — a7! a(n — 12°22! an) 

Boo = (G29 — aera, *ay2)—! = aoe! + arog eter" (1, — 2" 2201) ~* 412" 22 





Eqs. (9) will be employed to find Y*, and it is seen 
that a method must be devised for evaluating ay~' 
and a2~'. Now ay can be considered the sum of pi: 


Clearly Eq. (10) may be applied to a2~' and other 
reciprocations with obvious change of symbols. To 
illustrate the method of reciprocation and, in particular, 


and a tensor J, where J has n? tensor elements, the the scheme of partitioning to be used, the calculation of ' 

(n,n)th of which is K, and the others are nuil tensors. a~! is given in detail. Let Y; be expressed in par- 

The reciprocal of p; = Yi contains known tensor  titioned form: ( 

elements, and it is required to find a,,—' in terms of - ¥,.. ¥.. | 
—1 — -1 1= oc tan wa=- 

these and J. Thus an (p:1 + J). Now perform 7... 


the operations: 


pian = pi(pi + J) 
anpi~' = (a. + Do = I + Joi") 


where J is the unit tensor. 


an = ay (I + Joi") 


According to the subscript scheme, Y; has been parti- 
tioned between the (nm — 1) and mth compound indices, 
It is clear that Yi,, is a compound square tensor of 
(n — 1) ? elements and that Vig, and Ving are column 
and row tensors, respectively. By partitioning the 
tensor J between the same indices as Yi, the product 
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C—' may be found similarly from Eq. (10). 
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Omitting further intermediate steps, the admittance 
tensor, p*~! = Y* of the combined system can now be 
expressed. The compound tensor elements are ex- 
pressed in the simpler alternative forms of Eqs. (9) and 
are partitioned according to the scheme: 
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where 

D= (I+ (K’ — Liaunn"L) Vontinty)— (20) 
E= I+ (K — Lanatyinty Li) Vinal = (21) 


For convenience, the elements of a2.~!, omitted in the 
discussion, together with a repetition of a;,~', are given 
below. 


@1aa' = YViaa — Yian(Z + K Yinn)~'K Yina (22a) 
Qian? = Vian(I + K¥Yinn)~ (22b) 
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Onn) = Vina + K Yin)" (22d) 
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Examination of Eqs. (19) through (23) discloses 
that calculation of the Y* elements is not so formidable 
as it might first appear. For example, to compute a 
typical element, Y*12.»,, of Y*i2q, knowledge of K, 
K’, and L and the four original admittance tensors, 
Vinny Youtn(aty, Vian, and YVont+yo is sufficient. 
Table 2 summarizes the requirements. 
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Missing Elements are Supplied by Conditions of 
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APPLICATION OF THE THEORY 


Two structures, in which the degree of complexity 
prevents determination of their parameters, are to be 
connected in a known way. It is desired to predict 
the vibration response of various elements in the com- 
bined structure because of action of known forces 
imposed at selected points. Development of the theory 
shows that experimental determination of the original 
admittances provides sufficient information to solve 
the problem and that no intrinsic structural details are 
needed. . 

The treatment supposed that interconnection was 
made between one body of each structure, but this 
was a restriction employed for convenience only. 
Actually, a number of interconnections among bodies 
may be made and the resulting system analyzed, 
although the computations are considerably more 
involved. This more general case will be discussed 
presently, but now applications in which various simpli- 
fications demonstrate the flexibility and simplicity of 
the method can be considered. 


Vibration Isolation of Power-Plant Installations—General 


In aircraft, where the mass of the power plant is a 
substantial part of the airplane mass and where the 
air frame certainly does not act as a rigid body, con- 
ventional methods of calculating the desirable sus- 
pension flexibilities often lead to ineffective isolation 
of the power-plant exciting forces. Usually, the vibra- 
tion engineer must accept the air-frame and power- 
plant dynamic properties as they stand and restrict 
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his attention to means of connecting the two com- 
ponents. The method of analysis presented herein 
permits an indication of the complete airplane vibration 
characteristics; for, by applying calibrated excitation 
to the engine mount structure and determining the 
response of various fuselage components, the air-frame 
admittances, in the range of engine-propeller excitation 
frequencies, may be found. These, together with the 
power-plant admittances and selected values of sus- 
pension stiffness, enable the fuselage vibration to be 
predicted for known power-plant exciting forces. A 
final value of suspension stiffness may be chosen on 
the basis of the predicted overall performance. It 
should be noted that the equivalent viscous effect of 
structural and still air damping is implicitly taken 
into account; advanced instrumentation allows ac- 
curate determination of the amplitudes and phases 
of the admittances. 


Some Features of Vibration Isolation 


(a) In most applications, interest is centered on 
components of the structure to which the exciting 
system is appended. Accordingly, if S; is the original 
structure, the only significant part of Y* is the area 
Y* 2 of Eq. (19). 

(b) Frequently, the appended system consists of a 
few known elements, in which case the quan- 
tities Vorn4nin+1) and YVorn41, may be computed 
with greater facility than by experimental determi- 
nation. 

(c) Cases where the added system is merely a rigid 
body prove especially easy. Only Y*teq,+1 and 
Y*ea(nt+1) are required, which in turn demand experi- 
mental values of Vien and Vin»; Yon+1(nt+1) iS simply 
the diagonal tensor —(1/w?)|m;,)—'. Applications of 
this type are the installation of heavy machinery in 
buildings and ships and installation of aircraft engines, 
neglecting the effect of propeller-blade flexibility. If 
the propeller be considered as a rigid inertia, the effect 
of its rotation on the calculated engine admittances 
can be found. However, since the effect of spin is to 
couple the yaw and pitch coordinates, this calculation 
becomes rather lengthy. 

(d) An interesting application of power-plant vibra- 
tion isolation is one in which the vibration characteris- 
tics of the nonrotating propeller are included. The 
admittance tensor at the engine center of gravity may 
be found from experimental values of propeller ad- 
mittances and calculated engine admittances, some- 
what as discussed above. This total power-plant 
admittance is then used in conjunction with the air- 
frame admittances to provide the final result. It 
should be mentioned that the effect of rotation on the 
honrotating propeller admittance has as yet not been 
derived. 

(e) Often, specification of the system does not re- 
quire 6th order compound coordinates because of 


structural symmetry. When only two or three indi- 
vidual coordinates are needed to define each body, the 
simplification of calculation procedures is apparent. 

(f) A further reduction of effort accrues when the 
forces are simple or are compounded of only two or 
three simple forces. These conditions are by far more 
common in engineering applications than the general 
case. 


Additional Considerations 


It has been mentioned that the interconnection 
between two subsystems was restricted to the chain 
type in the interest of clarity. The method may be 
extended to include a seemingly hopeless combination 
formed by connecting r bodies of S; with s bodies of 
So. 

Assign identifying subscripts to the coordinates in 
such sequence that the impedance tensor, p*, may be 
written 
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Clearly, a similar analysis can now be applied with 
further partitioning made between the indices (m — r) 
and (n —r +1), and (m+ s)and(m+s+1). Eggs. 
(19) through (23) then give the components of Y*, 
where the subtensors Y*,,,, and Y*2(n+1(n+1) Change 
from unit order to r and s and the others are suitably 
revised. 

The schematic solution may be written Y* = 





a r s 6 
[12,... a-riliren... 2 Une... rslinsed).. inepl 









(a-r) 
(a-rei) 









(25) 





PRX A 
LNWS 


A study of athe vibration characteristics of multi- 
engined aircraft installations is particularly susceptible 
to analysis by this more general method. Ordinarily 
the air-frame system would be combined with the first 
power plant by means of Eq. (19) to produce the first 
modification. This process would demand reoperation 
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with Eq. (19) on the first modification and the second REFERENCES 
power plant, and so on, until all power plants had been + von Karman, Th., and Biot, M. A., Mathematical Methods in 


“installed.” Eq. (25), however, gives the complete Engineering, 1st Ed., pp. 370-378; McGraw-Hill Book Company, 
Inc., 1940. 


airplane admittance in one operation. The calcula- “sit we A., Coupled Oscillations of Aircraft Engine-Propeller 
tous involved are practicable — the interconnec- systems, Journal of the Aeronautical Sciences, Vol. 7, No. 9, pp. 
tion tensors K, K’, L,,, and L,, are diagonal. 376-382, July, 1940. 


. . . 3 i f } 4 . 7- e = 
Evaluation of the interconnection tensors for both f any Sineny Senet eens ay Meteniin, gp. TS; Tae 
Wiley & Sons, Inc., 1939. 


chain and general connections may usually be facilitated ‘ Frazer, R. A., Duncan, W. J., and Collar, A. R., Elementary 
by a change to relative coordinates and use of the trans- Matrices, pp. 112-113; Cambridge University Press, London, 
formation tensor.’ The inconvenience of operating 1938 

with complex admittance tensors may be avoided by mee uae = en pp. oe 

replacing them by tensors wth twice their num- 7 Concardia, c, The Use a Tensors in Mechanical Engineering 
ber of rows and columns containing only real quan- pyosjems, General Electric Revue, July, 1936. 
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At a meeting initiated by Dr. J. C. Hunsaker and Dr. Th. von Karman, acting as joint secre- 
taries of the Fifth International Congress for Applied Mechanics, a committee was formed for the wh 
organization of the Sixth Congress. In line with the decision reached at Cambridge in 1938, it is ' 
proposed that the Sixth International Congress for Applied Mechanics be held in Paris, from Sep- gay 
tember 22 to September 29, 1946. qu 
The invitations to the Congress are extended on behalf of 1’Académie des Sciences de 1|’Institut fre 
de France, la Direction des Relations culturelles, le Centre national de la Recherche scientifique, at 
l'Institut de Méchanique de la Faculté des Sciences de Paris, la Société francaise des Mécaniciens, of 
and l’Association technique Maritime et Aéronautique. - 
The Congress wilf meet at the Sorbonne. 7 
| The Congress will be divided into the following Sections: sho 
| I. Structures. Elasticity. Plasticity. -* 
II. Hydro- and Aerodynamics. Hydraulics. a 
III. Solid Dynamics. Vibration and Sound. Friction and Lubrication. . 
IV. Thermodynamics.: Heat Transfer. Combustion. Fundamentals of Nuclear Energy. ing 
Besides the papers presented in these Sections, a number of General Lectures will be given on stre 
| subjects of current interest. The titles of these Lectures will be made known in a later notice. wak 
| Those who desire to become members of the Congress are requested to inform the Secretary ra 
the 


General as soon as possible of their intention to attend. Communications are to be addressed to the 
Secretary General of the Sixth International Congress for Applied Mechanics, Institut Henri- 
Poincaré, 11, rue Pierre-Curie, PARIS (V°). 

The Proceedings of the Congress, containing the minutes of the sessions, the texts of the Gen- 
eral Lectures, and the papers presented in the Sections, will be published with the least possible de- T 
lay. The members of the Congress will have the right to buy these volumes by subscription. The 


surplus copies will be sold later on at a higher price. is © 
The payment of an inscription fee of 300 francs is required from the members of the Congress; bee 

this fee can be paid immediately or not later than on the opening day of the Congress. Ladies — 
accompanying members will not be required to pay such fee, except those who wish to address the R 
meetings or to submit papers. * 
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Vortex Frequency and Flow Pattern in the 
Wake of Two Parallel Cylinders at Varied 


Spacing Normal to an Air Stream’ 


HERMAN M. SPIVACK 
National Bureau of Standards 


ABSTRACT 


The frequency of vortex formations in the flow about a pair of 
parallel cylinders, at various separations, is investigated with the 
hot wire anemometer. A résumé of the literature on vortex fre- 
quencies is given with a view to relate the studies of obstacles of 
diverse shapes to the present case. The law of dynamical simili- 
tude is developed for individual solids, for the efflux through a 
rectangular orifice, for secondary eddies, and for the two-cylinder 
system. 

When two cylinders are separated by a gap just smaller than 
the diameter, instability occurs. At larger gaps the cylinders 
behave like independent bodies. As the gap is decreased, the 
main sequence of frequencies (probably generated on the outer 
sides of the cylinders) changes from a value corresponding to a 
single cylinder, of diameter d, to a value associated with a solid 
body of breadth equal 2d (at zero gap). A minimum frequency 
is reached at G/d = 1/1, where G is the gap distance. For gaps 
smaller than the critical, a formula is given to represent the main 
sequence: 


fd/U = 0.09 + {[(G/d) — */s]/\(G/d) + 2)}* 


where f is the frequency and U is the velocity. 

The picture is complicated by vortexes generated within the 
gap. At spacings less than one-half diameter, a low gap fre- 
quency is found. Between one-half and one diameter a high gap 
frequency is present, decreasing to the independent cylinder value 
at the critical spacing G/d = 1. In addition, doubled frequencies 
of the main sequence are also found, generally on the centerline 
of the wake, attributed to overlapping vortexes from the external 
sides of the cylinders. 

For all values of the gap, the dimensionless parameter fd/U is 
shown to be independent of Reynolds numbers greater 
than 15,000. Below this value, fd/U rises at first and 
then drops sharply, the maximum apparently depending on 
the gap. 

The velocity distribution in the wake was obtained by travers- 
ing with a pitot static tube. In a plane four diameters down- 
stream from the center of the system, it was found that the 
wakes of both cylinders are mixed for gaps below the critical and 
separate for larger gaps. At G/d less than about one-half, 
the central pressure drop is great enough to cause the flow to 
reverse. 


INTRODUCTION 


= FORMATION OF VORTEXES in the flow down- 
stream from obstacles immersed in moving fluids 
is of interest in fluid dynamics and acoustics and has 
been the subject of numerous investigations. 
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* A thesis on this topic has been submitted to the Graduate 
School of the Catholic University of America in partial fulfillment 
of the requirements for the degree of Master of Science. 


Early Investigations 


Mallock,' in 1907, was apparently the first to explain 
that the aeolian tones investigated by Strouhal, Ray- 
leigh, and others were the result of vortexes being shed 
alternately in the wake of cylindrical wires subjected to 
an air current. In 1908, Bénard? published the results 
of an extended experimental investigation of alternating 
vortexes. 


Stability of the Vortex Street 


In a remarkable investigation, von Kérmaén* has 
derived the criterion for the stability of the two rows of 
vortexes forming a two-dimensional street and has also 
indicated the relationship between the resistance of 
immersed objects and the vortexes associated with the 
flow. The stable type of vortex formation is an asym- 
metrical pattern that is referred to in the literature as 
the von K4rm4n or Bénard vortexes. In this case the 
vortexes are detached alternately from each side of an . 
obstacle, and the centers of the successive eddies main- 
tain a constant separation. 


Under ideal conditions, if the distance between the 
centers of the vortex rows is and the vortexes in each 
row are spaced a distance /, the spacing ratio h// has 
the constant value 0.281 for stability. In the experi- 
ments of von K4rm4n and Rubach‘ made on cylinders 
in water, photographic methods gave the values 0.305 
and 0.280. Tyler obtained 0.308 for cylinders. In the 
Fage and Johansen experiments® with a flat plate in the 
wind tunnel, this ratio is found to change from 0.250 at 
a downstream distance of five times the plate width, to 
0.52 at a downstream distance of 20 times the width of 
the plate (the channel width was 14 times the plate 
width). Change in the spacing ratio is accompanied 
by a diffusion of vorticity, and it is problematical that 
vortexes, as such, can be said to exist far downstream. 
The actual arrangement of the vortexes can be consider- 
ably affected by external conditions, such as the proxim- 
ity of the walls of a channel, the viscosity of the 
fluid, and the turbulence of the stream. Rosenhead® 
has shown that the presence of channel walls causes h// 
to increase downstream. An analysis by Hooker’ 
attributes the downstream variation in the spacing 
ratio principally to the viscosity of the fluid, while 
Dryden and Heald* have found that changing the 
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turbulence in the stream modifies the vortex configura- 
tion. 


Frequency of Vortex Formation 


The frequency at which two-dimensional cylindrical 
vortexes are shed is a linear function of the velocity 
over the entire range of the stable formation except at 
Reynolds numbers below 1,000. Otherwise, it has been 
demonstrated” that the linear range holds up to 
Reynolds numbers of 5 X 10°. Expressing the fre- 
quency in the form of the dimensionless parameter 
fd/U, for a circular cylinder, von Karman‘ obtained 
0.207 and 0.198; Relf and Simmons,’ 0.18 to 0.20; 
Blenk, Fuchs, and Liebers,'' 0.207. Tyler’? found 
somewhat lower values. In the plane of symmetry be- 
hind an immersed body it is possible to measure the 
effect of vortexes from each row overlapping to produce 
a double frequency. 


Edge Tones 

The presence of alternating vortexes may be noted 
for an entirely different flow condition. This was 
shown, for example, by Benton" in the study of edge 
tones produced when a thin plate is attacked by a jet 
of air. The jet issuing from a narrow rectangular aper- 
ture enters a body of still air with no contraction. The 
vortexes generated yield for the ratio #// values ranging 
from 0.35 to 0.276, depending upon the width of the slit. 
These values are established indirectly and their deter- 
mination involves acoustic considerations. Earlier 
investigators have examined various acoustical aspects 
of this problem. Carriere'* has obtained interesting 
photographs of such vortexes. 


Secondary Vortexes 


A brief reference to the so-called secondary vortexes 
as distinguished from alternating vortexes may be de- 
sirable. These are sometimes produced at the borders 
of wakes of semiplanes (i.e., a sharp-edged plate ex- 
tending into the fluid from one of the boundaries) 
placed in a current. Such secondary vortexes have 
also been observed with cylinders," in which case they 
begin at a Reynolds number of about 1,250 and also 
tend to reach a limiting upper frequency and may be- 
come coincident with the principal vortex formation. 
They are of relatively small dimensions and high fre- 
quency. The most extended study of this subject was 
undertaken by Crausse, whose results are included in an 
informative discussion by Camichel and Escande."* 
Photographs have also been taken by Prandtl” of the 
formation of secondary vortexes. 


The Present Problem 


While numerous investigations have been conducted 
upon individual obstacles having continuous surfaces 
and diverse shapes, few studies have been made of the 
vortex pattern associated with more complex forms. 
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The extension of the study of alternating vortexes may 
naturally lead to multiple cylinders, the simplest ar- 
rangement being the case of a pair of cylinders separated 
by various distances. This case, however, is not the 
simplest for a theoretical study. Various theoretical 
problems may find their guidance in the actual observa. 
tion of vortexes. In this paper the principal study will 
be of the frequencies appearing in the wake of two 
cylinders, dependence of the frequency on velocity, and 
the effect of the separation of the cylinders. 


SYMBOLS 


= total cross-stream breadth of an obstacle 

= diameter of a cylinder (!"/s in.) 

frequency, cycles per sec. 

= gap between, or separation of, the adjacent walls of the 
cylinders 

= Reynolds number (Ud/v) 

= Strouhal number (fd/U) 

“u = local velocity in the x-direction at any point in the flow, 

ft. per sec. 

Au = fluctuations of the mean velocity 

U = free stream velocity 

= coordinate direction parallel to the wind stream 

coordinate direction normal to the wind stream 

= kinematic viscosity 


Q» ac 
ll 


nw 


vee 
Il 


LAWS OF DYNAMICAL SIMILITUDE FOR VORTEX 
FREQUENCIES 


Similarity Law for Cylinders 


It was observed by Strouhal’* that the frequency of 
the sound emitted by a wire of diameter d, moving in air 
at a velocity U, obeys, within limitations, the relation 


S = fd/U (1) 


where S is constant. Bénard’ referred to the quantity 
Sas the Strouhal number. Lord Rayleigh’ examined 
the data of Strouhal and, assuming negligible compressi- 
bility of the air, derived the expression for alternating 
vortexes: 

S = 0.195{1 — (20.1/R)] (2) 


where R is the Reynolds number Ud/v. This is the 
law of dynamical similitude for the vortexes leaving the 
cylinder and is in good agreement with the formula 
obtained in an investigation by Tyler:'” 


S = 0.198{1 — (19.7/R)] (3) 


These formulas are valid at low Reynolds number down 
to the point where the eddies that are formed do not 
detach from the cylinder, so that the Karman trail no 
longer exists. Thom” finds this limit to be about R = 
30, and Camichel'® gives the value R = 48. The 
physical significance of this criterion is that for flow at 
lower Reynolds.numbers the eddies that are formed do 
not detach from the cylinder, and any disturbances are 
damped out. 

The Strouhal number has been shown to be constant 
for the range R = 10° to R = 5 X 10°. This is approx- 
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imately the range in which the drag of a cylinder has a 
constant value. Different investigations are in fair 
agreement regarding the numerical magnitude of S for 
a cylinder. Dryden and Heald*® obtained 0.19; 
Camichel!® and his collaborators, 0.20; Relf and 
Simmons,’ from 0.18 to 0.20. Measurements made 
with a single cylinder in connection with the present 
investigation yielded 0.207, which is the value given by 
Blenk, Fuchs, and Liebers.!! 


Similarity Law for Prismatic Solids 


For prismatic solids the dynamical similitude law is, 
of course, similarly expressed as for the cylinder. 
Taking 5 as the breadth of an object whose length tra- 
verses a two-dimensional flow, the Strouhal number is 
given by fb/U and the Reynolds number by )U/v. 
There also appears to be a range of R for prismatic 
solids wherein S has a constant value. Camichel’® 
gives the following: for a rectangular body with 
depth three times the breadth, S = 0.156; for a square 
rod, S = 0.130; for a rectangular block terminating in 
an equilateral wedge, S = 0.235; for a wedge tapering 
downstream, with the base three times the depth, 
0.255. 

With flat plates fb/U is nearly constant for angles of 
incidence exceeding 30°. Fage and Johansen® obtained 
an average value of 0.148; von Karman,’ 0.144; 
Tyler,'* 0.158; and Blenk, Fuchs, and Liebers,'! 0.180. 
For airfoils, Tyler found the value 0.15, whereas Blenk, 
Fuchs, and Liebers found 0.21. For the shapes of a 
wedge, an ogive, and an extended ogive, respectively, 
the values 0.238, 0.271, and 0.254 are given by Fage and 
Johansen,*' who are of the opinion that the frequency 
and spacing ratio depend upon the lateral spacing of 
the vortex sheets leaving the body and, to a lesser ex- 
tent, upon the form of the body. ‘The relation of Sto R 
shows marked difference for solids of various shapes, 
but the values obtained for similar shapes by different 
investigators show the order of dispersion to be almost 
as marked. 


Similarity Law for Jets 


Under certain conditions, vortexes are formed in the 
efflux from an orifice. For the vortexes associated with 
the flow in jets, dynamical similitude for the vortex 
frequency may be expressed in the form: 


fd/U = F,\(Ud/y) (4) 


where F; is a functional symbol, d is the width of the jet, 
and U is the velocity of efflux. In the experiments of 
Benton,'* when the edge plate is close to the slit, no 
tone is obtained. However, when the edge is moved to 
a greater distance ad, a note suddenly appears. For 
the determination of a, Benton gives the dynamical 
similitude law: 


(ao/d)? = 1.50 + (Br/ Ud) (5) 


where B is a constant approximately 2,000. In view of 


Konig’s formula 
= U/2ao (6) 


where f is the frequency of edge tone, Benton’s equation 
may be modified to yield a result conforming with Eq. 
(4). Putting S = fd/U, Konig’s formula becomes 


S = d/2ay (7) 


and, when this is introduced in Benton’s equation, Eq. 
(5), the new form is 


1/45? = 1.50 + (Bv/Ud) (8) 


This relationship may be depended upon to give satis- 
factory results for small Reynolds number. The curve 
giving the relation between S and Reynolds number 
may be obtained from the data in Fig. 5 of Benton’s 
paper. This curve as derived from the original data is 
given in Fig. 2. The value of constant S obtained by 
this indirect means is approximately 0.32, but it must 
be remembered that the data are for low Reynolds 
numbers. 


Similarity Law for Secondary Vortexes 


The secondary vortexes have been best observed 
when they arise from a semiplane. The law of dy- 
namical similarity for the frequency of the secondary 
vortexes takes on a different form in this case. Con- 
sider a point at a distance X from the edge of a semi- 
plane and located within the same plane as the latter. 
The velocity at X will be U, and, since the dimensions ° 
of the boundaries may be said to have no effect on the 
flow, the Reynolds number will be referred to the 
distance X. The similarity law will then be given in 
the form: 


U/V fv = Fo(XU/v) (9) 


Camichel’s’ analysis of the experiments of Crausse 
and Baubiac gives (at X = 6) approximately 


U/W fv = 28 (10) 


For the secondary vortexes in the wake of a cylinder the 
appropriate law should be of the form 


U/V/ fv = F(dU/v) (11) 


where d is the diameter of the cylinder, 


Dynamical Similarity for Two Cylinders 


Depicting the Flow.—The vortexes formed behind two 
parallel cylinders may assume a configuration that is 
not readily susceptible of an analysis. In the absence 
of precise information regarding the structure of the 
wake, we may conceive of a particular arrangement of 
the vortexes. The data obtained by observations may 
then be examined in the light of such conception. One 
vortex street may be set up at the external sides of the 
two cylinders corresponding to a flow about a solid body 
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with the same overall dimensions. Another type of 
street may be set up about each cylinder, corresponding 
to the conventional flow about an individual cylinder. 
A third type of vortex formation corresponding to flow 
through a narrow rectangular orifice may also be en- 
visioned. Complication of the flow patterns may be 
caused by interaction of the various vortex streets and 
by secondary vortexes. For simplicity of treatment, 
we may restrict ourselves to just two types of vortex 
flows, internal and external. The frequency of the 
alternating vortexes generated at the external borders 
of the two cylinder system must now be expressed by 


fd/U = F,(Ud/v, G/d) (12) 


where F, is a functional symbol. The diameter, d, of 
each cylinder, and G, the gap distance between them, 
are the only geometrical quantities governing the flow. 
When G exceeds a certain critical value to be considered 
later, the frequencies become practically independent 
of the gap. The frequency of the vortexes generated 
at the internal boundaries of the cylinders—that is, at 
the gap—are likewise expressed by a form 


fa/U = F,(Ud/v, G/d) (13) 


Here again, when G exceeds the same critical value 
mentioned above, the frequency becomes independent 
of G. For separations of the cylinders greater than the 
critical gap, Eqs. (12) and (13) become identical. As 
the size of the gap is changed, marked differences must 
occur in flow through the gap. Such changes may be 
discontinuous and reflect on the manner in which the 
vortexes are formed. Visual or photographic observa- 
tion could serve to resolve this problem, but in the 
absence of such information we may further contem- 
plate the details of the internal flow. 


Internal Vortex Picture 


When the gap is zero the system behaves as a single 
solid of breadth equal 2d. Since there is no flow 
through the internal boundaries, the vortexes are 
generated only on the external sides. It may be possi- 
ble, however, for a stationary vortex pair to be set up 
in the downstream cavity formed by the cylinders, and, 
when the outer vortexes are detached, the stationary 
pair may oscillate or be dislodged. Proceeding to the 
case of a small gap, a stable jet may be issued with a 
parabolic velocity distribution. Being laminar, the 
jet may widen downstream and become bifurcated upon 
encountering the principal flow in the wake. The 
separating branches of the jet would be endowed with 
vorticity of opposing sign, forming two cylindrical 
eddies. As the gap is further increased in size, the jet 
may issue in a cone of uniform velocity between two 
vortex sheets which may unfold into an asymmetrical 
vortex street. In this state the mechanism of forma- 
tion of the inner and outer vortexes differs but little. 
The difference, if it exists, would be in the effective 
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Fic. 1. Traversing mechanism. 

circulation of the inner and the outer eddies. As a 
result of the conditions of flow for various gap values, 
the law given in Eq. (13) may show various discontinui- 
ties. The secondary vortexes may, of course, be present 
in any of the cases discussed above, but, since they may 
become confluent with the more prominent vortexes and 
thus indistinguishable, no further statement can be 
made at present. 


EXPERIMENTAL SETUP 


The study was conducted in the 4'/2-ft. wind tunnel 
of the National Bureau of Standards. The cylinders 
used were of commercial brass tubing 1°/, in. in diam- 
eter, 1/,.-in. wall thickness, mounted to extend com- 
pletely across the 4'/2-ft. working section of the wind 
tunnel. The setting of the gap between the cylinders 
was held to within 0.04 in. in the plane of the measure- 
ments and to closer tolerance for the smaller spacings.* 
The cylinders were rigidly clamped to steel angles at the 
walls of the tunnel and did not vibrate perceptibly for 
any of the speeds encountered. The flow in the empty 
working section is of low turbulence (see footnote on 
page 299) and an examination with threads disclosed it 


* A slight bend in one of the cylinders (0.04 in. in 41/, ft.) was 
neglected as having no measurable effect, since random inter- 
change of the cylinders throughout the investigation did not give 
rise to any attributable discrepancies. 
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TABLE 1 
Frequencies of Vortexes Downstream from Two Parallel Cylinders. Velocity = 40 Ft. per Sec. 
Main Sequence 

Gap, ——_——— Frequencies ——-——~ Doubled Gap Other 
In. G/d Coupled Uncoupled Frequencies Frequencies Frequencies 

0 0 35.4 71 

1/32 0.023 36.1 72 

1/16 0.045 36 70 

3/32 0.068 33 66 

1/8 0.091 30.6 

2/8 0.182 31.7 15.5 
3/8 0.273 34 68 17 20-22 
4/8 0.364 34 

5/8 0.455 35.7 70 16,114 16, 22 
6/8 0.545 39.5 

7/8 0.636 42 80 115 

8/8 0.727 46.6 90 110-116 70, 80, 90 
9/8 0.818 50 100 100 
10/8 0.909 5&8 74 16, 25, 84, 90 
11/8 1.000 73-77 150 36, 93 
12/8 1.091 2 
13/8 1.182 72 142 
33/16 1.500 74 
22/8 2.000 71 140 
33/8 3.000 73.5 142 
44/8 4.000 7s 146 
55/8 5.000 73 145 

72 144 


Single cylinder 72 


to be as uniform as visual observation will reveal. An 
absence of frequencies capable of interfering with the 
cylinder measurements, in the flow in the empty tunnel, 
was also indicated by an exploration with the hot wire. 

Measurements were made of the velocity pressure 
distribution across the wake at several stations down- 
stream from the cylinder system. To accomplish this, 
the flow was traversed with a pitot static tube, using the 
traversing apparatus shown in Fig. 1, especially con- 
structed for the purpose. The accuracy of such meas- 
urements was within 2 per cent. Frequency measure- 
ments were made in similar explorations of the flow 
with the National Bureau of Standards Hot Wire 
Anemometer,*? with which the velocity distribution was 
also checked in a few cases. Readings were taken in a 
plane parallel to the windstream approximately midway 
between the ends of the cylinders. Thus, wall interfer- 
eice and end effects were reduced to a minimum. The 
hot wire mount was kept to small dimensions and the 
carriage of the traversing apparatus was, 1n every case, 
held at no less than 15 diameters downstream from the 
point at which the readings were being taken in order to 
minimize interference between the carriage and the 
cylinders. Vibration of the hot wire mount, Fig. 1, 
was not easily eliminated because of the narrow depth 
cantilever, unsupported for about 20 in. However, 
vibration was not detected in the direction normal to 
the two-dimensional plane of flow, and it is in this plane 
that the wire is most sensitive. An oscillation did 
appear at a very low frequency (about 3 cycles) at a 
wind speed of 40 ft. per sec., and, since it occurred in 
the direction parallel to the longitudinal axes of the hot 
wire and of the cylinder system, it is believed to be of 
minor importance. 

To obtain the frequency with which a hot wire 
(platinum, 0.0003 in. in diameter, 0.2 in. long, heated 


by a current from 50 to 100 milliamp.) was periodically 
cooled when subjected to an eddying flow, the output of 
the hot wire anemometer was impressed on the hori- 
zontal plates of an oscilloscope. The output of a 
variable frequency oscillator (Hewlett Packard) was 
impressed on the vertical plates and tuned until the 
screen of the oscilloscope showed a simple closed figure. 
The preferred figure is a circle, representing a 1:1 
correspondence of the frequency picked up in the. 
stream to the calibrating frequency. By tuning to 
Lissajou figures of multiples, a check on the frequency 
is obtained. Fage and Johansen® indicate that rela- 
tively steady readings could be obtained in the region 
outside the wake. However, with practice, the observer 
was able to recognize frequencies present in all but the 
most turbulent sections of the wake. Extraneous fre- 
quencies generated in line voltages, amplifying circuits, 
and other pickup could be reduced to negligible propor- 
tions by reducing the oscilloscope amplification and in- 
creasing the heating current to the hot wire. Thus the 
hot wire pickup was generally made overwhelmingly 
predominant on the screen. The error of the frequency 
measurements was within 2 per cent at all speeds, ex- 
cept in the critical region (see “Results of Frequency 
Measurements’ below) where it may be double this 
value. 


RESULTS OF FREQUENCY MEASUREMENTS 


Classification of the Frequencies 


In exploring the flow behind and between the cylin- 
ders, a number of different frequencies were obtained 
for the same air speed and cylinder spacing. Usually 
these frequencies occurred at different positions in the 
field of flow, but different frequencies were also found 
simultaneously in the same place. Figs. 6 to 15 show 
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the positions at which the measurements were made for 
the various gaps and the precise distribution of the 
frequencies found at an air speed of 40 ft. per sec. 

The observed values are collected in Table 1. These 
are classified as follows: The frequencies measured in 
the general domain of the cylinder wakes are referred to 
as the “main sequence.’’ When the flow about one 
cylinder has an effect on the flow about the other, as at 
small separations, the main sequence will be said to 
consist of “‘coupled frequencies.”” When the cylinders 
are separated by a larger gap, where the flow about one 
does not appear to influence the frequency at the other, 
the main sequence will be called “uncoupled fre- 
quencies.” 

Another group of frequencies, these measured 
distinctly in the vicinity of the gap, are listed 
under “gap frequencies.” Miscellaneous frequencies 
that could not readily be classified were also recorded, 
and these are listed as ‘‘other frequencies.’’ Corre- 
sponding to the main sequence is a group that registered 
frequency values twice as high as the latter; these are 
given as ‘‘doubled frequencies.” 


Doubled Frequencies Explained 


The existence of the doubled frequencies can be ex- 
plained as follows: When the vortexes detached from a 
cylinder in a von K4érman train have developed to the 
magnitude at which they reach or overlap the center- 
line of flow, each side of the street will contribute equally 
to the measurement at that point and effect a fre- 
quency twice as high as that obtained in the general 
flow (as given in the “‘Introduction’’). Reference 5 also 
reports such observations. In the case of the cylinders 
touching, the double frequency is, in fact, found on the 
centerline and can be explained by the overlap of the 
outer sides of the vortex street. This is also true for 
gaps of !/3: and */3 in., as in Figs. 6 and 7. For gaps of 
18/4, 33/5, 44/5, and °*/, in. (Figs. 10, 13, 14, and 15), the 
doubled frequencies occur only at the centerlines of the 
individual cylinders, as in the case of a single cylinder. 
At the intermediate gaps no simple statement can be 
made. Whether the doubled frequencies, which are in 
such cases more widely distributed over the field, owe 
their origin to an overlap of the vortex formation from 
the two cylinders, cannot be said with certainty. 


Main Sequence of Frequencies 


A graphic representation of the data of Table 1 is 
given in Fig. 3. Disregarding the curve corresponding 
to the doubled frequencies, it is seen that the observa- 
tions taken at values of G less than '"/; in. distribute 
themselves in three distinct groups. All the fre- 
quencies at gaps larger than ''/, in. appear to fall on a 
line of constant ordinate, indicating that the flow about 
the cylinders is uncoupled, and each cylinder generates 
vortexes of like nature. For the ordinate having the 
value f = 72 cycles per sec., the cylinder diameter ''/. 
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in., and the velocity of air 40 ft. per sec., the correspond- 
ing Strouhal number is fd/U = 0.207. 

If the assumption is made that the main sequence is 
due to vortexes formed in an external street (and this 
assumption is not unfounded), a satisfactory explana- 
tion of the phenomena may be arrived at. When the 
cylinders are touching, there is only one dominant fre- 
quency and it must be external. For this case 


f(2d/U) = 2f4/U = 0.21 


and the frequency of the vortexes is about 36 cycles per 
sec. As the gap is enlarged, the frequency of the ex- 
ternal vortexes changes from the initial value of 36 toa 
final value of approximately 72. The flow modification 
at the external borders is a continuous process up to a 
gap of '/, in., where an unstable distribution seems to 
produce several prominent frequencies. Thus, the 
lower branch of the two curves that join at G = ""/s in. 
represents the variation of the external frequencies as a 
function of gap distance. This is the main sequence 
curve for the coupled frequencies and exhibits a mini- 
mum at G = '/sin., rising in an approximate parabola as 
it approaches G = "/, in. in one direction and more 
sharply toward G = '/3 in. where the frequency is 
slightly higher than at G = 0. 


Vortexes at Small Gaps 


The decrease in frequency from G = 0 to G = !/; in. 
may be explained by considering the flow through the 
gap of closely spaced cylinders as of a nature to delay 
the separation of the vortexes at the external borders. 
Since the boundary layer separates from the surface of a 
cylinder in the form of a vortex sheet, which at Rey- 
nolds number encountered in this investigation is in 
the form of a discreet vortex, a definite quantity of 
fluid must be available at each air speed to maintain 
this periodic detachment. If a leak of the fluid in the 
boundary layer occurs (i.e., at very small spacing) the 
process of vortex formation at the external sides may be 
retarded. Prandtl” has shown that removal of fluid 
from the boundary layer can delay or prevent separa- 
tion. 


Vortexes in the Gap 


Returning to the remaining groups of frequencies, 
the highest branch, descending as it approaches the 
critical value, was tabulated as a gap frequency be- 
cause it was recorded with some intensity in the gap. 
A peculiarity of this descending branch is that the fre- 
quencies are also noted in the general wake of the two 
cylinders. This could mean that the vortexes from an 
internal jet produce periodic pressure variations that 
cause velocity fluctuations. Another possibility sug- 
gested by the scattering of the readings is that these 
vortexes are unstable and meander under the action of 
the external pressures. That these internal vortexes 
are not quite similar to the edge-tone vortexes dis 








ob 
flo 
ga) 
otl 
ob 
ide 


Ot. 


Fis 
ing 


we 





‘espond- 


uence is 
and this 
explana- 
hen the 
ant fre- 


cles per 
the ex- 
36 toa 
ification 
up toa 
eems to 
us, the 

11/, in. 
eS asa 
quence 
a mini- 
abola as 
d more 
ency is 


= 1/, in. 
ugh the 
o delay 
orders. 
ace of a 
it Rey- 
nm is in 
tity of 
aintain 
| in the 
ng) the 
may be 
rf fluid 
separa- 


1eTICIes, 
ies the 
icy be- 
le gap. 
he fre- 
he two 
rom an 
is that 
y sug- 
t these 
tion of 
ortexes 
es dis- 





VORTEX FREQUENCY 295 








d=/88mm 
148 
LiF 
98 
79 


36 
20 





Ud 








290 +90 i'd fp ___iqoo__iapo__itpo _/f@ __‘/800 ___ 2090 
























Fic. 2. Vortexes in the jet from a rectangular orifice. 


cussed above results from a consideration of gap dimen- 
sions and measured frequency. A connection is sug- 
gested, however, if the constant value of 2.S, shown in 
Fig. 2, is considered to be too high for the orifice that is 
formed by two adjacent cylinders in an airstream. For 
a gap of °/,in., for example, Fig. 2 yields a frequency of 
about 140, whereas 113 was observed, but the dis- 
crepancy becomes progressively greater at smaller gaps. 
This example, however, assumes a linear extrapolation 
for Reynolds numbers from 2,000 to 28,000. 

Since it is believed that extraneous vibrations, as of 
the cylinders and apparatus, do not enter into the pic- 
ture, the distinctness and separation of the two branches 
of the gap frequencies suggest two possibilities. First, 
there may be two modes of vortex formation at the 
borders of the internal jet, or, second, as the gap dis- 
tance is increased two distinct régimes are established, 
depending upon the cylinder separation. As may be 
seen from Fig. 3, none of the high gap frequencies were 
obtained below G = °/s in., which indicates that the 
flow does not support these frequencies at the smaller 
gaps or that they are too weak to be detected. On the 
other hand, none of the low gap frequencies were 
observed at gaps larger than °/s in., which supports the 
idea of a discontinuity at that point. 


Other Frequencies 


The low gap frequencies indicated by the squares of 
Fig. 3 may be contrasted with a number of other read- 
ings, of the same order of magnitude, which could not 
be obtained distinctly.. These readings, tabulated as 
“other frequencies,’’ are shown by x's in Fig. 3. They 
were not sharply defined but spread over a band of 
several cycles and were found rather close to the cylinder 
at gaps above, as well as below, °/sin. Since such read- 
ings were also obtained with a single cylinder (and were 
absent with both cylinders removed), this may be the 
limiting value of the secondary vortexes mentioned by 
Courregelongue.'® The lack of such data for cylinders, 
however, prevents the making of conclusive state- 
ments, 
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diameter cylinders). 


DEPENDENCE OF FREQUENCY ON REYNOLDS NUMBER 


The pertinent data of Fig. 3 are reproduced in Fig. 4 
in the form of dimensionless parameters. This repre- 
sentation is for a Reynolds number of 28 X 10°, taking 
the kinematic viscosity of air as 16 X 10~° ft.* per sec. 
It is desirable to examine the data to determine whether 
or not there is any observable effect when the Reynolds 
number of the flow is varied. In particular, it is im- 
portant to know whether the critical value of G/d de- 
pends upon the Reynolds number. To determine this 
point the values of the Strouhal number for values of the 
Reynolds number from 5,000 to 93,000 were superim- 
posed on Fig. 4. The result was that all the values 
above R about 15,000 fell within a few per cent of each 
other and were indistinguishable on the curve and, 
hence, omitted. 

At R = 12,000 the curve was about 10 per cent 
higher, and at R = 5,000 the curve was about 10 
per cent lower than the mass of measurements at 
greater values of R, for all gaps. This indicated that a 
scale effect exists for the lower velocities. 
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Fic. 6. Frequency distribution about two cylinders, touching 
and separated by !/32 in. 


Linear Behavior of Strouhal Number 


To examine this condition further, Fig. 5 has been 
prepared, where the data relate principally to the main 
sequence of external vortex frequencies. As indicated 
in the “‘Introduction,’’ the limits of the stable vortex 
formation is also the region of constant fd/U for a given 
body. 





3/32 and 3/, in. apart. 


Frequency distribution about two cylinders spaced 


For this condition to be fulfilled the frequency 
must be a linear function of the velocity of the airstream. 
Such behavior of the frequency is shown for all cylinder 
spacings in Fig. 5, where, for R > 15,000, the lines 
drawn through the observed points are practically 


parallel to the ordinate axis. 
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Fic. 8. Frequency distribution about two cylinders spaced 
8/s and °/s in. apart. 


in. for the Reynolds numbers below 15,000, the Strouhal 
number rises and then drops sharply toward the low 
end. While sharp definition of the lowest frequencies is 
difficult, it nevertheless appears that a transition in the 
flow has occurred. This agrees with the discussion of 
the preceding paragraph. The departure occurs within 
the usually stable range for a single cylinder. In this 
respect, it may be deduced that the two-cylinder system 
is an anomalous case. 


An Empirical Formula for the Frequency 


When the two cylinders are far removed from each 
other (uncoupled), one will produce a vortex street as if 
the other were absent and the Strouhal number may be 
expressed by fd/U for either. On the other hand, as the 
gap approaches zero, the limiting frequency of the total 
system must approach a value for which S is given by 
f(2d/U). For the intermediate case, between these 
two limiting conditions of the external vortex forma- 
tion, a law of the form 


fa/U = F,(G/d) (14) 


may be used to express the frequency parameter. This 
isin accord with Eq. (12) above, since S has been shown 


Fic. 9. Frequency distribution about two cylinders spaced’ 
11/, and '*/s in. apart. 


to be largely independent of Ud/v. To express this 
lower branch of the main sequence curve of Fig. 4, the 
following empirical formula is given: 


fa/U = 0.09 + {((G/d) — '/s]/[(G/d) + 2]}? (15) 


Eq. (15) is valid for gap parameters less than the critical 
value, and the accuracy with which it approximates the 
observations is shown in Table 2. 


VELOcITY MEASUREMENTS 


Velocity pressure measurements were made for the 
several cylinder spacings with the pitot ‘static tube 
shown in Fig. 1. Readings were taken in a plane nor- 
mal to the flow mainly four diameters downstream from 
the center of the system. In this plane (x/d = 4), the 
pitot static tube was generally out of the region of 
severe pressure changes, stagnation, or reversal of flow 
behind the cylinders. The dynamic pressures were 
obtainable, therefore, without error caused by direc- 
tional flow around the pressure tubes. 

Figs. 11 to 15 show the results of such measurements. 
In Fig. 16 u/U, is plotted for several spacings both 
smaller and larger than the critical value of G/d = '/,. 
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13/, and *3/s in. apart. 


At x/d = 4 the wakes behind the cylinders are definitely 
separate for gaps greater than the critical. At smaller 
gaps the wakes have definitely merged to exhibit pres- 
sure distributions like those at a greater distance down- 
stream from a single continuous obstacle. Such 
behavior is significant in supporting the frequency 
observations above and below the critical gap 
value. 

At G/d = ‘**/,, the velocity along the central axis of 
flow has reached the free stream velocity in the plane 
x = 4d. It is natural, therefore, to suppose that the 
flow from each cylinder is here, independent of the 
other. At G/d = */, the central velocity does not 
quite reach the free stream value so that slight coupling 
is assumed. The central pressure has dropped still 
farther at G/d = '*/\, although mixing of the wakes has 
not yet occurred. At G/d = °/,,; a single wake pattern 
exists with minimum pressure at the axis of symmetry. 
While single wake patterns are expected below G/d = 
0/,,, the central pressure region is flat at °/;,; but be- 
comes steep immediately at smaller gaps until, as at 
4/,,, the flow has reversed direction and cannot be 
measured precisely with the pitot static tube. Mix- 
ing of the wakes from each cylinder occurs earlier 


























cylinders spaced °/s; in. apart. 


for smaller spacings, but for measurements made 
directly behind the system individual wakes were de- 
tected at all but the smallest gaps. At G/d = */, the 
wakes were found to mix between one and two diameters 
downstream. The positions shown in Figs. 6 to 15 are 
the approximate limits beyond which mixing or ex- 
cessively weak signal prohibited reading of the fre- 
quencies. The outer boundaries of the wake can be 
discerned for all gaps by pressure measurements, at 
least as far downstream as x/d = &, although mixing 
has obviously obliterated the initial vortex frequencies 
(see Fig. 99, reference 23). 

In the wake of a single cylinder, reversal of flow is 
characterized by stationary vortexes attached to the 
rear of the body.* Ata single cylinder, the dead space 
(the region of reversal of flow) is wedge-shaped, with 
the base of the wedge at the cylinder and the vertex 
about 2d downstream. For two cylinders close to- 
gether (G < d) but not touching, there are two such 
dead space wedges, even for a gap '/id, one for each 
cylinder. But the vertexes are about ld downstream, 
so that—e.g., for x = '/sd—the width of each dead 
space is smaller than in the case of the single cyl 
inder. 
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Fic. 12. Velocity and frequency distributions about two 


cylinders spaced '°/s; in. apart. 


At the downstream half of each cylinder up to a dis- 
tance of 0.25 in., the hot wire exploration shows a high 
degree of turbulence, which seemed to be more notice- 
able than for single cylinders. 

Measurements of the mean speeds were obtained in 
traverse of the wake with the hot wire anemometer. 
Resulting readings were substantially as obtained with 
the pitot static tube for a few spacings observed, as in 
Figs. 6 to 15. As applied, the single hot wire is not 
sensitive to directional changes in the x-y plane so that 
only the magnitude of the mean speed and not the 
velocity is disclosed. When the single hot wire is used 
in coordination with the pitot static tube, reversal of 
flow, as well as the magnitude, may be detected. The 
circled values (in the upper portion of Fig. 11) of the 
speed measured with the hot wire are shown to be in a 
negative direction by the pitot static tube measure- 
ments. 


Vorticity 


The N.B.S. Hot Wire Anemometer circuit is equipped 
to indicate the value (squared) of the fluctuations about 
the mean speed. This arrangement has been used? to 
determine the intensity of the turbulent fluctuations in 
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Fic. 13. Velocity and frequency distributions about two 


cylinders spaced 2*/3 in. apart. 


the airstream. While quantitative results can be ob- 
tained by the procedure outlined in reference 24, Table 
IX, the circuit is admirably suited to furnish a direct 
qualitative indication of the vorticity distribution in the 
wake. The magnitude of the fluctuations is a function 
of the vortex strength at the point that readings are 
being taken, and the square of the velocity fluctua- 
tions has the dimensions of vortex strength per unit 
area. 

In a von Karman vortex street, the fluctuations have 
definite frequencies. In the free stream the existing 
vorticity appears as “isotropic turbulence,’’*> in which 
discreet vortex frequencies cannot be detected but a 
more or less continuous spectrum of eddy frequencies of 
much smaller scale and intensity than the von Karman 
vortexes exists. This requires the wake of upstream 
obstacles to be thoroughly mixed and diffused before 
reaching the point in question.* (For details see 
reference 26.) Upstream from the two cylinders or far 
to one side—i.e., in the free stream—the von Karméa 
vorticity is not of measurable proportions. It becomes 


‘ 


* In the N.B.S. 4!/2 ft. wind tunnel working section the turbu- 
lence at 40 ft. per sec. is of the order of 100(u’? + v’? + w’*)/3U = 
0.02 per cent. 
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Fic. 15. Velocity and frequency distributions about two 


cylinders spaced *°/s in. apart. 


measurable as the wake is approached. The fluctua- 
tions increase as the wake is entered, rising violently at 
the center of the wake whether the latter is formed by 
two individual streets or by a single coupled pattern. 
This is illustrated in the upper portions of Figs. 11, 12, 
13, 14, and 15, where (Az)? is plotted to an arbitrary 
scale dependent on the dimensions of the wire, the 
heating current, and the velocity. 
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Distribution of velocities in wake of two cylinders 
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CONCLUSIONS 


For single obstacles immersed in a flowing fluid, it 
may be said, in general, that the Strouhal number 
fb/U is a constant, where b is the cross-stream breadth 
and independent of the Reynolds number for the turbu- 
lent flow régime. ‘The same statement may be applied 
to the case of two cylinders, except that the constant 
value fd/U changes for (G/d) <1. For Reynolds num- 
bers, Ud/v, less than 15,000, a nonlinear relation be- 
tween frequency and velocity indicates a transition in 
the state of flow. G/d = 1 is the critical spacing above 
which the cylinders behave independently. At smaller 
values of the gap, individual von Karman vortex 
streets are not detached from each cylinder; instead, a 
coupled vortex street is formed by the overall system 
along with certain internal vortexes that depend upon 
the gap. The internal vortexes are probably generated 
in vortex sheets bordering a jet that emanates from the 
orifice formed by the adjacent sides of the two cyl- 
inders. 

A discontinuity may be said to exist in the internal 
vortexes at G/d = °/, below which a low frequency, 
and above which a high frequency, is found. The fre- 
quencies of the coupled external vortex street are given 
approximately by 


fa/U = 0.09 + {[(G/d) — ¥/s)/[(G/d) + 2)}" 


This formula holds up to G/d = '/,, where the fre- 
quench exhibits an instability. A minimum frequency 
at G/d = '/y is explained by delay of separation of 
vortexes, caused by a condition analogous to a leak in 
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the boundary layer at a solid body, when the two 
cylinders are close together. 
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ABSTRACT 


If the possibility of a variation of load factor with altitude is 
admitted, then the cabin portion of a pressurized aircraft is sub- 
jected to two sets of loads, both of which vary with altitude. 
The first set of loads is composed of the loads due to the flight load 
factors; the second set, of all the loads imposed by internal pres- 
sure in the cabin. 

In the section of this paper titled ‘‘Variation of Load Factor 
with Altitude” it is shown that the flight load factors will probably 
decrease with altitude. Since pressure loads will increase with 
altitude, it follows that various structural components in the 
cabin region will be critically loaded at what may be designated as 
the critical structural altitude (i.e., since the structural members 
are loaded by one set of loads decreasing with altitude and another 
increasing with altitude, there must exist some altitude at which 
the state of stress in a particular member becomes critical. It is, 
of course, possible that this altitude will be sea level or possibly 
the aircraft absolute ceiling). In general, the critical altitude 
will be different for each structural component. 

Of the conclusions reached in this paper, there is one of particu- 
lar interest because it represents a departure from accepted 
practice. This conclusion is: 

In general, it is not possible to assume that diametral pressure 
loads are carried by hoop tension in the skin unless a careful in- 
vestigation proves that there is no other structure capable of 
carrying these loads (this holds for fuselages of any cross section). 


VARIATION OF LOAD FACTOR WITH ALTITUDE 


HERE ARE TWO TYPES of flight loadings to which 

an aircraft structure may be subjected. These are 
(1) maneuvering loads and (2) gust loads. To deter- 
mine a true picture of the variation with altitude of 
either of these loading types is beyond the scope of this 
paper, and unnecessary as far as the following discus- 
sions are concerned. The only fact that it is desired to 
establish is that the magnitude of both these loading 
types probably will decrease in some fashion with 
altitude. 


Maneuvering Loads 


Since maneuvering load factors are in general arbi- 
trarily assigned, it is difficult to state exactly how they 
will vary with altitude. It is known, however, that the 
higher up an airplane flies, the less maneuverable it be- 


comes. This follows from the fact that the maneuver- 
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ability of an aircraft is largely dependent on the effec- 
tiveness of its control surfaces which, in turn, is de- 
pendent in part on the air density (which decreases with 
altitude). To make this point clearer, it is only neces- 
sary to imagine a projectile flying in some manner 
through space with only conventional tail surfaces 
available for control. It is obvious, since the air 
density would be zero, that the projectile would be un- 
controllable (that is, unmaneuverable). 
On the basis of the above, 


tm = f(p) = fh) 
where 


Nm = maneuver load factor 
f(e) = function of air density 
f(h) = function of altitude 


Since there are many other factors that influence the 
value of the maneuver load factor, it should be recog- 
nized that the only point brought out above is the de- 
pendence of the maneuver load factor on altitude. 

It now remains to demonstrate that the gust load 
factor also varies in some fashion with altitude. 


Gust Loads 


In discussing gust loads for pressurized cabin air- 
craft, it is necessary to keep in mind that the primary 
purpose of pressurization is to permit the airplane to 
avoid turbulent air conditions, such as those encoun- 
tered in thunderstorms, fronts, and the like (i.e., over 
the weather flight). 

It is known that extremely severe gusts are encoun- 
tered even at the base of the stratosphere in thunder- 
storm cloud formations. However, it is not logical to 
assume that pressurized cabin aircraft will fly through 
this type of disturbance when these aircraft are specifi- 
cally designed to avoid this type of flight. 

Working on the sensible assumption that pressurized 
aircraft will avoid turbulent air conditions at all times 
except during take-offs and landings, and recalling that, 
except during storm conditions, gusts are associated 
with ground discontinuities (mountains, lake, shore, 
etc.), it is readily seen that the magnitude of gusts which 
a pressurized airplane may be expected to encounter 
will decrease in some fashion with altitude and become 
negligible in the stratosphere.! 

On the basis of the above, it appears that the gust 
velocity that should be used to design pressure cabin 
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aircraft is a function of altitude and decreases with in- 
creased altitude. 
The accepted formula? for gust load factor increment 


is: 
An, = (p/2)KUV(dC,/da)/(W/S) 

where 

An, = gust load factor increment 

p = air density 

K = a factor depending on wing loading 

U = gust velocity 

V = true air speed 

dC,/da = lift slope 

W/S = wing loading 


In the above formula, p, U, and V may vary with alti- 
tude. It follows, therefore, that An, = f(); further- 
more, An, is apparently a decreasing function of h, 
which eventually becomes zero. 


An Arbitrary Altitude Load Factor Reduction Formula 


In order to permit numerical results to be presented in 
the next section, the following arbitrary formula for 
load factor variation with altitude is assumed. 

nr, = 1 + B Ang (1) 
where Amp = load factor increment at sea level and 

0.6942h 7 
— ——_; 0 < h < 35,382" 

35,332 (2)* 
8B = 0.3058; h > 35,332’ 





B=1 


TERMINOLOGY 


P, = internal pressure running load (lbs. per in.) 

3,¢ = angles 

Xm, Xp, X» = unknown bending moment, axial load, and shear 
load, respectively 


p = pressure 

Po = pressure at sea level 

x, y, = rectangular coordinates 

a.b = sides of plate 

t = thickness 

v = Poisson’s ratio 

u,v,w = components of displacements 

W = running load at section 

E = Young’s modulus 

oj,» = normal stresses due to inertia or flight load 


Tn, Tp, 01 = Stresses due to pressure 
1, Ft,» Ft,» 01, = total stress 


= VEI/P, 





] 

I = moment of inertia 

P. = axial load in beam column or frame 
L = beam length 

M = moment 

A, Ay = area stringer, area frame cross section 
r = frame radius 

P = concentrated inertia load on frame 


* Up to the base of the stratosphere, 8 as given by this equation 
is merely a linear approximation of the density ratio p/p); above 
35,332 ft. 8 was assumed to be constant. 


DETERMINATION OF CRITICAL STRUCTURAL ALTITUDE 
AND DISCUSSION OF STRESS ANALYSIS PROBLEMS 
PECULIAR TO THE DESIGN OF PRESSURIZED CABINS 


Hoop Tension and Membrane Tension Stresses in Skin 


Although in earlier works* on the topic of this paper, 
hoop tension in the skin was considered as the major 
factor resisting radial pressure due to pressurization, it 
is not possible to tacitly assume that this is always the 
case. Consider the following argument pertaining to a 
fuselage of circular section. 

(a) For hoop tension to exist in the skin the skin 
must be continuous. 

(b) If hoop tension in the skin is to be of importance, 
internal framework must have less resistance to radial 
and circumferential, deformation than the skin. 

(c) Since the skin of the modern metal monocoque 
fuselage is composed of numerous small panels riveted 
to each other and to the internal framework and, since 
in the cabin or pressurized portion of the fuselage the 
skin is of necessity cut by windows, doorways, and 
escape hatches, the skin may not be considered per- 
fectly continuous. As a ‘consequence, condition (a) 
above is violated to some degree. It follows that the 
skin may be incapable of carrying hoop tension loads. 

(d) In the cabin portion of a fuselage large concen- 
trated loads are applied to the structure (wing reactions, 
landing gear reactions, seat loads, cargo loads, and the 
like). These loads require the use of strong, rigid bulk- 
heads and heavy stringers. Additional heavy frames 
are required in many instances to prevent buckling of 
the skin since the greatest fuselage bending moments 
occur at the cabin portion. 

Before assuming that the skin (in hoop tension), 
rather than the internal framework, carries the internal 
pressure loads, an elastic analysis should be made which 
proves that the skin is less elastic than its supporting 
structure. 

Although the skin may not be capable of carrying 
diametral pressure loads in hoop tension, it will in any 
event serve as a means of transmitting the pressure 
loads to the internal structure. Since the skin thick- 
ness will generally be small in comparison to the edge 
dimensions of the panel considered, membrane theory 
may be used in the analysis of the skin. 

If one assumes that the boundary conditions provided 
by the riveted attachments of the skin to the internal 
framework may be approximated by considering the 
skin panel to be simply supported along the edges but 
restrained at the edges from moving in the plane of the 
skin, then the development presented on pages 345 and 
346 of reference 4 may be extended to panels of rec- 
tangular shape.t 

The results of such an extension for a Poisson’s ratio 
of vy = 0.3 (Timoshenko chose 0.25) for aluminum are 

t It is interesting to note that there is an arithmetical error in 
Timoshenko’s work on page 346.‘ In Eq. (204), for wy» the 
numerical constant should be 0.817 rather than 0.802. 
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TABLE 1* 
orzaty = b/2,x =0,a; cyatx =a/2,y = 0,bfory = 0.3 
or and a, vs. a/b 





a/b o7/(a*p*E/t*) oy/(a*p*E/t*) 
1/s 0.575 0.258 
1/s 0.585 0.339 
2/s 0.560 0.400 
3/4 0.542 0.427 
1 0.460 0.460 





* Curves of o;/(a*p?E/t?) and o,/(a*p*E/t?) are plotted in Fig. 
3 


summarized in Table 1 and presented graphically in 
Fig. 1. 

The stresses o, and a,, as given in the table and in the 
figure, are the stresses at the midpoints of the panel 
sides, since the calculations made indicated that the 
stresses attained their maximum values at these points. 

In applying either Table 1 or Fig. 1, it should be noted 
that the origin of coordinates is at a corner of the plate 
whereas in reference 4 the origin was taken as the center 
of the plate; hence the dimension a of this paper is 
twice as large as a in reference 4. 


Effects of Pressurization on the Ability of a Skin Panel to 

Carry Shear or Compression 

In the event that a panel buckles under either shear 
or compression, the buckled pattern is generally in the 
form of a number of waves that require both inward and 
outward displacement of the panel. (This is also 
true in the case where the pattern involves more than 
one panel.) Internal pressure retards the formation of 
the inward buckles and aids the formation of outward 
buckles. On this basis it might be assumed that 
internal pressure would leave the buckling load un- 
altered, since the amount of additional energy required 
to produce the inward buckles appears to be offset by 
the smaller amount of energy required to produce out- 
ward buckles. 

A theoretical investigation of the problem is beyond 
the scope of this paper; however, the limited amount of 
test data available to date seems to indicate that the 
buckling load is not reduced by internal pressure. 

On the basis of the above and the fact that the dis- 
cussions in ‘‘Variation of Load Factor with Altitude”’ 
indicate that compression and shear loadings will 
probably decrease with increased altitude (and, conse- 
quently, with increased internal pressure), it may be 
stated that the critical altitude for panels subject to 
shear and compression will probably be sea level. It 
should be noted when considering compressive stresses 
in the skin produced by the flight loads that the mem- 
. brane tension stresses due to pressurization will be re- 
lieving stresses. 

Extending the above discussion, it may be concluded 
that internal pressure has the net effect of increasing 
the effective width of the skin under compression. 
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Fic. 1. Membrane tension stresses at the midpoints of the 


panel sides. The direction of the stress is shown in the small 
inset figure. 


Effects of Pressurization on the Ability of a Skin Panel to 

Carry Tensile Loads 

‘The foregoing discussion deals with two of the three 
types of stresses to which the skin is subject in perform- 
ing its structural functions (i.e., shear and compression 
stresses). In the discussion of these two types of stress, 
the conclusion was reached that the critical structural 
altitude for a panel loaded by either type was probably 
at sea level. This statement was permissible, since 
loads causing compression or shear decrease in propor- 
tion to the load factor (which should decrease with alti- 
tude) and therefore are greatest at sea level. It will be 
shown below that this statement does not apply to skin 
panels under tensile loads. 

Tensile loads in the skin of an airplane fuselage are 
imposed, in general, by external or inertia loads that 
vary with the load factors, and therefore with altitude. 
The skin of a pressurized cabin is loaded also by mem- 
brane tension stresses, and longitudinal stresses caused 
by its action as a carry-through between the fore and 
aft pressure-sealing bulkheads. As a consequence, the 
skin is loaded by one set of loads increasing with altitude 
and another set decreasing with altitude. ‘The critical 
structural altitude of a typical skin panel is determined 
below. 

On the basis of Eq. (1) the stresses due to inertia or 
external loads are given as follows. 


o, = K,(1 + Ang) (3) 


where K, = stress at sea level/(1 + An). 
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From Table 1, the stresses due to normal pressure are 


on = Ky(p)” (4) 
where K, for a square plate would be 
K, = 0.46(a2E/t?)'” (5) 
For longitudinal stresses: 
o, = Ki(p) (6) 
where K, = the cross-sectional area of the cabin at a 


point considered divided by the total area of the skin 
and the stringers at the cross section and o, = total 
stress = 0; + on + CO}. 


o, = K,(1 + AnB) + K,(p)" + Kip) (7) 


To simplify calculations it is assumed that p (the in- 
ternal pressure) increases linearly with altitude, be- 
tween h = O and h = 35,332 ft. The error introduced 
by this assumption is in the order of 25 per cent at 
17,000 ft., but this error, as is shown later, does not 
affect the calculation of the critical structural alti- 
tude. 

The linear variation of p is expressed below (assuming 
sea-level pressure is maintained in the cabin). 


p = 11.289h/35,332; 0 < h < 35,332 ~— (8) 
In terms of Eq. (8), Eq. (7) may be written 


0.6942h = 4 ‘ 
35,332 


11.280h \"” 11.280h 
- a 35,332 ) +s ( 35,332 ) 
For o, maximum, 0o,/O0h = 0; whence 


= ( 311K, Y’ 0 
1.962 AnK,; — 31.95K,/ ’ 


From examination of Eq. (9), it is seen that if 
(1.962 AnK, — 31.95K,) is negative, then the critical 
altitude is at sea level. Furthermore, if this equation 
yields a value of h greater than 35,332, which is beyond 
the domain of applicability of the equation, it is appar- 
ent that the critical altitude will be the absolute ceiling 
of the aircraft. j 

For purposes of illustration the critical altitude / 
will be calculated for a typical skin panel. Values for 
An, K,;, K,, and K,, are determined below. 

Take An = 3.061. - K; is assumed to be 8,000 Ibs. 
per sq.in. 

To determine K,, assume that a = 10 in., E = 10° 
lbs. per sq.in., 4 = 0.064 in. ‘Then for a square panel 


10? X 10° 
0.064? 


= Ki(1+ An — 


< h S 35,332 (9) 





/s 
K, = 0.46 ( ) = 2,875 (Ibs. per sq.in.)”” 
For the determination of K,, a circular fuselage 12 ft. 
in diameter is assumed having a skin thickness of 0.064 
in. The area of the skin and stringers is assumed equal 
to twice the area of the skin. 


K, = [w(72)*]/["2(144)(0.064)] = 281.5 


Substituting the above values of An, K,, K;, and K, 
into Eq. (9), it is found that h = 12,200 ft. 

To evaluate the error introduced by the assumption 
of a linear pressure variation, the variation of total 
stress (o,) with altitude is determined from Eq. (7) using 
the true pressure variation, and the results plotted in 
Fig. 2. From examination of this figure it is seen that 
the critical structural altitude is 12,200 ft., which 
checks the results obtained from Eq. (9). 

It will be noted from examination of Fig. 2 that there 
is little variation in o, between oj,;,- ANG Cima. This is 
the result of the choice of K;, K,, and K,, and does not 
indicate that the variation of o, with altitude is small in 
all cases. 

In the calculation of ¢, vs. 4, no ultimate factor ot 
safety was assumed, since a true picture of stress varia- 
tion was desired. At this point it is of interest to note 
that when dealing with a structure possessing a non- 
linear load-stress relationship, it is desirable to base all 
margins of safety on the ratio of allowable yield stress 
to limit stress rather than the ratio of allowable ulti- 
mate stress to ultimate stress. 
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Fic. 2. Variation of total tension stress in skin panel with 
altitude. Note that the maximum stress intensity occurs at 
12,200 ft. 


Effects of Internal Pressure on the Load-Carrying Ability 
of the Longitudinal Stiffeners 
In general, the functions of longitudinal stiffeners re- 
sult in the stiffeners being in a state of either pure ten- 
sion or pure compression, if no internal pressure exists. 
In addition to the pure tensile or compressive stresses 
that normally exist, bending stresses due to lateral 
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loads are induced in the stiffeners when buckling of the 
skin, bulkheads, or of the stringers themselves occurs. 

To resist these lateral buckling loads and also to re- 
sist tendencies toward premature general instability 
failure, stiffeners are designed to have a large amount of 
flexural rigidity in comparison to the skin or plating. 
As a consequence of this, stringers function as beams 
under lateral load. 

Lateral loads on stiffeners or stringers are generally 
ignored in analysis; however, in designing stringers for 
a pressurized cabin, it is impossible to ignore lateral 
loading, since lateral loads are built up in direct propor- 
tion to the internal pressure. The effects of lateral 
pressure loading of the stringers due to internal pressure 
will now be considered. 

Distribution of Pressure Load on Stringer 
facilitate computations, the following 
assumptions are made regarding pressure load distribu- 
tion on the stringers. 

(1) The total normal plate load is divided equally 
between the four boundaries: i.e., if the panels on both 
sides of the stringers are the same size, then each panel 
causes a load abp/4 + abp/4 = abp/2, where a and b 
are the plate edge dimensions. 

(2) The load distribution is in the form of a sym- 
metric triangle. 

(3) Corner forces are negligible. 

(4) The stringers are simply supported and may be 
considered to be composed of separate parts of length 
equal to the bulkhead spacing. 

Stiffener Loaded by Axial Tension—In the event that 
the axial load is tension, it is conservative to ignore the 
secondary moment. 

Ignoring the secondary moment 


M = (WLx/4) — (Wx?/3L); 0 < x < L/2 
Maz, at x = L/2 = (WL?/12) 


In order to 
simplifying 


where 
L beam length 
W = running load at apex of triangle 
M = bending moment 


o, = bending stress = (WL?/12/)c 
where 


distance from stringer neutral axis to stringer 
edge 
I = stringer moment of inertia 


c 


In addition to the bending load due to normal pres- 
sure, the stiffener is subject to the same axial tension 
loads as the skin, namely: 

Co; = K,(1 oe AnB) and Co; = K(p) 


hence 


=H +o+o,= i + Ki(l + Ang) + Ki(p) 


From “Distribution of Pressure Load on Stringer,”’ 


WL/2 = abp/2 or W = abp/L; whence 
o, = K,(p) + Kil + AnB) + Kip) (10) 


where A, = abLc/12I/. 

A plot of o, vs. h (considering true pressure varia- 
tion) is presented in Fig. 3 using the following values for 
the constants involved. 


K; = 8,000 Ibs. per sq.in. 





An = 3.061 
K, = 281.5 
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Fic. 3. Variation of maximum tensile and compressive stresses 
in stringer and effective width of skin (compression) with alti- 
tude. 


Fig. 3 indicates that the maximum stringer stress 
occurs at sea level. This result should not be inter- 
preted to mean that all stringers subject to tension are 
critical at sea level, inasmuch as a different choice of 
constants could lead to a critical altitude either at the 
aircraft ceiling or, in special cases, at some intermediate 
altitude. 

Stiffener Loaded by Axial Compression—It can be 
shown that the maximum moment for a simply sup- 
ported beam under axial load and subjected to a sym- 
metric triangular loading occurs at the center of the 
beam and has the value 


M = —(2Wj7*/L)(tan L/2j) + Wj? (11) 
where W and L are defined above and 


J? = EI/P, 
i axial load 
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Fic. 4. Variation of maximum tensile stress in the circular bulk- 
head (shown at the top of the figure) with altitude. 


The axial load P, acting in a stringer under compression 
is the difference of the compression stress o; and the ten- 
sion stress ¢;, multiplied by the area of the stringer and 
its associated effective width of skin, i.e.: 


fe = (o; eS, aA, 
where A, = area of stringer and effective width of skin. 
Substituting the equations for o, and @;: 
P, - [K,(1 + Ang) wis K(p))A; 


The curve of o,, vs. h is plotted in Fig. 3 for the values of 
K;, An, K,, A;, W, L, I/c, a, and 6 assumed below: 


K, = 8,000 Ibs. per sq.in. 
An = 3.061 

K, = 281.5 

A, = 0.60 sq.in. 

lf ="'/, 

c =1 

W = abp/L = 10p 

f = El/P, = 10°/3P, 


From examination of Fig. 3, it is seen that the critical 
altitude is sea level as in the case of o,. As in the pre- 
vious case, the conclusion should not be drawn that the 
critical altitude will always be sea level. 


Effects of Internal Pressure on the Load-Carrying Ability 
of Bulkheads Loaded in the Plane of the Bulkhead 
The analysis presented in preceding sections has 
shown that a critical structural altitude (which may be 
sea level or any other altitude up to and including the 


aircraft’s absolute ceiling) exists for both stringers and 
skin. It is possible to generalize and state that a criti- 
cal altitude exists for any structural component sub- 
jected to the simultaneous action of a variable internal 
pressure and variable flight or inertia loads. The proof 
of this statement is of an intuitive nature and, conse- 
quently, the statement may be considered axiomatic. 

A method of determining the critical structural alti- 
tude of a typical bulkhead is now considered. The 
bulkhead to be analyzed 1s shown in Fig. 4. 

Using the method of analysis of reference 5 it can be 
shown that for the problem of Fig. 4: 


Fy ioe 
X» = p,r 
where 


X,’ = axial tension at cut due to internal pressure 


p, = running load due to internal pressure 

r = radius of ring (ring is circular) 

X,, = 0.1901Pir 

X, = 0.5000P; Pinternal loads at cut due to P; 


X, =0 


Making the assumption that the ring is critical in bend- 
ing and that shear stresses may be ignored, the critical 
stresses occur at y= 0 (see Fig. 4 for ¢). 

At this point, Mar. = 0.6901Pir and axial tension = 
py. The total stress is then: 


, = (0.6901 P,rc/T,) + (P,r/Ay) (12) 
where 
c = distance from the neutral axis to outermost 
fiber 
ZI; = moment of inertia of frame cross section 


A, = area of frame cross section 


Letting 0.6901rc/T, = Ki r/Ay = x, and v4 = KP, + 
Kypy. 

Since P; is an inertia load, its variation may be ex- 
pressed in terms of Eq. (1): 


o1 = K,Po(1 + AnB) + Kp, (13) 


Assuming a bulkhead spacing of 10 in., p, = 10p; where 
p is the differential pressure. 


Let 
K,Po = 8,000 Ibs. per sq.in. 
K, = 100/in. 
An = 3.061 
then 
o, = 8,000(1 + 3.0618) + 100), (14) 


To illustrate the variation of o,, with altitude, values 
of o,, have been calculated at various altitudes and the 
results plotted in Fig. 4. From examination of this 
figure, it is seen that o,, has a maximum value of 32,650 
Ibs. per sq.in. at 4,000 ft. 
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CONCLUSIONS 


(1) The load factor that an airplane may be expected 
to develop is apparently a function of altitude. How- 
ever, an equation expressing the exact variation cannot 
be derived without considerably more study. 

(2) Under the influence of a varying load factor there 
exists a critical structural altitude for each structural 
component of the pressurized portion of the airplane. 

(3) In general, it is not possible to tacitly assume 
that diametral pressure loads are resisted by hoop 
tension in the skin. 

(4) On the basis of available information it appears 
that the critical structural altitude for skin subject to 
shear or compression and the simultaneous action of 


internal pressure is at sea level. (Further research is 


indicated. ) 
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Calculation of Flow of Air and Diatomic 
Gases 


CHESTER W. SMITH* 


General Electric Company 


ABSTRACT 


Measurement of airflow is frequently required in connection 
with aeronautical projects. Conventional practice, as exempli- 
fied by the A.S.M.E. Test Code for Flow Measurement, provides 
two different formulas (or sets of formulas) for the calculation of 
flow from pressure and temperature measurements made with 
standard nozzles, orifices, and the like, the proper formula for use 
in any particular case depending upon whether the pressure ratio 
across the restriction is greater or less than the critical 
value. 

In most cases of aeronautical interest the formula for sub- 
critical (or restricted) flow is applicable. However, the exact 
formula for this case is quite cumbersome, and often requires 
much laborious calculation. Without sacrifice of exactness, the 
much simpler formula for critical flow may be generalized and 
used in all cases merely by including a ‘‘restriction factor” that is a 
function of pressure ratio only, and so may be tabulated once and 
for all. This paper proposes such a formula and presents a 
table of restriction factors applicable to air and other diatomic 


gases. 


SYMBOLS 


| 


= throat area 
discharge coefficient 

= specific heat at constant pressure 
specific heat at constant volume 
= acceleration of gravity 

= ratio of specific heats (= ¢,/c,) 
constants (Tables 2 and 3) 
restriction factor (Table 1) 

= absolute pressure 

Pi — be 

pressure ratio (= py,/p») 

specific weight 

gas constant 

absolute temperature 

specific volume 

V velocity 

w = gravimetric flow 

X = isentropic (or adiabatic) factor (= r —»/* — 1) 
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Subscripts 


0 = at selected standard or reference conditions 

at initial, or high-side, conditions 

at final, or low-side, conditions 

at discharge, or throat, conditions (i.e., just preceding 
throat area A) 

corresponding to critical pressure ratio 

based on total pressure, total temperature 


whom 
“ou i 


t 


(The subscript ¢ is added to a numerical subscript, as py, 714, 
%. Without the subscript ¢ it is to be understood that pressures 
and temperatures are static values and that specific volume and 
specific weight are based on static pressure and temperature.) 
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INTRODUCTION 


oe OF AIRFLOW is frequently necessary 
in connection with tests of aircraft power plants, 
accessories, and the ducts required in the airplane in- 
stallation. The usual method of making such measure- 
ments is by use of standard nozzles, orifices, venturi 
meters, or similar restrictions. From the measured 
pressures and temperatures the flow is calculated by 
well-known formulas that have been derived in many 
published papers and textbooks and which are summar- 
ized in reference 1. These formulas are also used to 
calculate the sizes of ducts and openings needed for 
specified flows and to determine the maximum flow 
capacity for given conditions. 

It is customary to distinguish sharply between 
different types of flow of a compressible fluid. If the 
total pressure and temperature just preceding the re- 
striction are constant, the flow per unit area gradually 
increases as the absolute pressure on the downstream 
side decreases, until a certain ‘‘critical’’ pressure ratio is 
reached at which it is a maximum. The flow at this 
pressure ratio is termed “‘critical.”” Further decrease 
of the downstream pressure does not change the flow, 
which under this condition is still termed “critical” or 
sometimes ‘“‘supercritical.”” When the downstream 
pressure is higher than that corresponding to the critical 
ratio, the flow is a function of this pressure and is termed 
“subcritical.” The words ‘“restricted’’ and ‘‘nonre- 
stricted,”’ or ‘‘affected’’ and “‘unaffected,’’ or ‘‘sub- 
sonic’ and “‘supersonic’’ are often used instead of ‘“‘sub- 
critical’’ and “‘critical,’’ respectively. 

The exact formula for subcritical flow is quite cum- 
bersome, and approximations have been proposed by 
different writers.» * However, for any group of gases 
having approximately the same ratio of specific heats or 
isentropic exponent, it is possible, by tabulation of a 
certain function called the “restriction factor’’ (perhaps 
“subcritical flow factor’’ would be a better name, but 
“restriction factor’ is the term originally applied and 
now generally adopted), to use without approximation 
what is essentially the simple formula for critical flow 
over the entire range of pressure ratios. Such a table 
(Table 1) is given here for air and diatomic gases, which 
have approximately the same ratio of specific heats, 
taken as 1.3947 for reasons explained later in more de- 
tail. A similar table could be easily made up for other 
gases if the need should arise. 
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TABLE 1 i Deron EY ae: 
Values of Restriction Factor N for Air and Diatomic Gases 
KCApyN A = throat area pu = high-side pressure (total) 
a C = discharge coefficient pf, = low-side pressure (static) 
VT iu K = constant (Table 2) Ty = high-side temperature (total) 
N = restriction factor w = gravimetric flow 


To obtain volumetric flow at any specified pressure and temperature, multiply w by the specific volume at that pressure and tem. 
perature. 

N is based on isentropic exponent 

k = 1.3947, or (k — 1)/k = 0.283 





Pu/ pe 0 1 2 3 4 5 6 7 8 9 
1.0000 0 0.006540 0.009249 0.01133 0.01308 0.01462 0.01602 0.01730 0.01850 0.01962 
1.0001 0.02068 0.02169 0.02265 0.02358 0.02447 0.02532 0.02616 0.02696 0.02774 0.02850 
1.0002 0.02924 0.02996 0.03067 0.03136 0.03203 0.03269 0.03334 0.03397 0.03460 0.03521 
1.0003 0.03581 0.03640 0.03698 0.03756 0.03812 0.03868 0.03922 0.03977 0.04030 0.04082 
1.0004 0.04134 0.04186 0.04236 0.04287 0.04336 0.04385 0.04433 0.04481 0.04529 0.04576 
1.0005 0.04622 0.04668 0.04713 0.04758 0.04803 0.04847 0.04891 0.04935 0.04978 0.05020 
1.0006 0.05063 0.05105 0.05146 0.05187 0.05228 0.05269 0.05309 0.05349 0.05389 0.05429 
1.0007 0.05468 0.05507 0.05545 0.05583 0.05621 0.05659 0.05697 0.05734 0.05771 0.05808 
1.0008 0.05845 0.05881 0.05917 0.05953 0.05989 0.06024 0.06059 0.06094 0.06129 0.06164 
1.0009 0.06198 0.06233 0.06267 0.06301 0.06334 0.06368 0.06401 0.06435 0.06468 0.06500 
1.0010 0.06533 0.06566 0.06598 0.06630 0. 06662 0.06694 0.06726 0.06757 0.06789 0.06820 
1.0011 0.06851 0.06882 0.06913 0.06944 0.06974 0.07005 0.07035 0.07065 0.07095 0.07125 
1.0012 6.07155 0.07185 0.07214 0.07244 0.07273 0.07302 0.07331 0.07360 0.07389 0.07418 
1.0013 0.07447 0.07475 0.07503 0.07532 0.07560 0.07588 0.07616 0.07644 0.07672 0.07699 
1.0014 0.07727 0.07754 0.07782 0.07809 0.07836 0.07863 0.07890 0.07917 0.07944 0.07971 
1.0015 0.07997 0.08024 0.08050 0.08077 0.08103 0.08129 0.08155 0.08181 0.08207 0.08233 
1.0016 0.08259 0.08284 0.08310 0.08335 0.08361 0. 08386 0.08412 0.08437 0.08462 0.08487 
1.0017 0.08512 0.08537 0.08562 0.08586 0.08611 0.08636 0. 08660 0. 08685 0.08709 0.08734 
1.0018 0.08758 0.08782 0.08806 0.08830 0.08854 0.08878 0.08902 0.08926 0.08950 0.08973 
1.0019 0.08997 0.09020 0.09044 0.09067 0.09091 0.09114 0.09137 0.09160 0.09184 0.09207 
1.0020 0.09230 0.09253 0.09275 0.09298 0.09321 0.09344 0.09366 0.09389 0.09412 0.09434 
1.0021 0.09457 0.09479 0.09501 0.09524 0.09546 0.09568 0.09590 0.09612 0.09634 0.09656 
1.0022 0.09678 0.09700 0.09722 0.09744 0.09765 0.09787 0.09809 0.09830 0.09852 0.09873 
1.0023 0.09895 0.09916 0.09937 0.09959 0.09980 0.1000 0.1002 0.1004 0.1006 0.1009 
1.0024 0.1011 0.1013 0.1015 0.1017 0.1019 0.1021 0.1023 0.1025 0.1027 0.1029 
1.0025 0.1031 0. 1033 0.1035 0.1038 0.1040 0.1042 0.1044 0.1046 0.1048 0.1050 
1.0026 0.1052 0.1054 0.1056 0.1058 0.1060 0.1062 0.1064 0. 1066 0.1068 0.1070 
1.0027 0.1072 0.1074 0.1076 0.1078 0.1079 0.1081 0.1083 0.1085 0.1087 0.1089 
1.0028 0.1091 0.1093 0.1095 0.1097 0.1099 0.1101 0.1103 0.1105 0.1107 0.1108 
1.0029 0.1110 0.1112 0.1114 0.1116 0.1118 0.1120 0.1122 0.1124 0.1125 0.1127 
1.003 0.1129 0.1148 0.1166 0.1184 0.1202 0.1219 0.1236 0.1253 0.1270 0. 1286 
1.004 0.1303 0.1319 0.1334 0.1350 0.1366 0.1381 0.1396 0.1411 0.1426 0.1440 
1.005 0.1455 0.1469 0. 1483 0.1497 0.1511 0.1525 0.1539 0.1552 0.1566 0.1579 
1.006 0.1592 0.1605 0.1618 0.1631 0.1644 0.1656 0.1669 0.1681 0.1693 0.1706 
1.007 0.1718 0.1730 0.1742 0.1754 0.1765 0.1777 0.1789 0.1800 0.1812 0. 1823 
1.008 0.1835 0.1846 0.1857 0.1868 G. 1879 0.1890 0.1901 0.1912 0.1922 0.1933 
1.009 0.1944 0.1954 0.1965 0.1975 0.1986 0.1996 0.2006 0.2017 0.2027 0.2037 
1.010 0.2047 0.2057 0. 2067 0.2077 0.2087 0.2096 0.2106 0.2116 0.2125 0.2135 
1.011 0.2145 0.2154 0.2164 0.2173 0.2182 0.2192 0.2201 0.2210 0.2219 0. 2228 
1.012 0.2238 0.2247 0.2256 0.2265 0.2274 0. 2283 0.2291 . 0.2300 0.2309 0.2318 
1.013 0.2327 0.2335 0.2344 0.2353 0.2361 0.2370 0.2378 0.2387 0.2395 0.2404 
1.014 0.2412 0.2420 0.2429 0.2437 0.2445 0.2453 0.2462 0.2470 0.2478 0.2486 
1.015 0.2494 0.2502 0.2510 0.2518 0.2526 0.2534 0.2542 0.2550 0.2558 0.2565 
1.016 0.2573 0.2581 0.2589 0.2596 0.2604 0.2612 0.2619 0.2627 0.2635 0.2642 
1.017 0.2650 0.2657 0.2665 0.2672 0.2680 0. 2687 0.2694 0.2702 0.2709 0.2716 
1.018 0.2724 0.2731 0.2738 0.2745 0.2753 0.2760 0.2767 0.2774 0.2781 0.2788 
1.019 0.2795 0.2803 0.2810 0.2817 0.2824 0.2831 0.2838 0.2844 0.2851 0. 2858 
1.020 0.2865 0.2872 0.2879 0. 2886 0.2893 0.2899 0.2906 0.2913 0.2920 0.2926 
1.021 0.2933 0.2940 0.2946 0.2953 0.2960 0.2966 0.2973 0.2979 0.2986 0.2992 
1.022 0.2999 0.3005 0.3012 0.3018 0.3025 0.3031 0.3038 0.3044 0.3051 0.3057 
1.023 0.3063 0.3070 0.3076 0.3082 0.3088 0.3095 0.3101 0.3107 0.3114 0.3120 
1.024 0.3126 0.3132 0.3138 0.3144 0.3151 0.3157 0.3163 0.3169 0.3175 0.3181 
1.025 0.3187 0.3193 0.3199 0.3205 0.3211 0.3217 0.3223 0.3229 0.3235 0.3241 
1.026 0.3247 0.3253 0.3259 0.3265 0.3271 0.3276 0.3282 0.3288 0.3294 0.3300 
1.027 0.3305 0.3311 0.3317 0.3323 0.3329 0.3334 0.3340 0.3346 0.3351 0.3357 
1.028 0.3363 0.3368 0.3374 0.3380 0.3385 0.3391 0.3397 0.3402 0.3408 0.3413 
1.029 0.3419 0.3424 0.3430 0.3435 0.3441 0.3446 0.3452 0.3457 0.3463 0.3468 
1.030 0.3474 0.3479 0.3485 0.3490 0.3495 0.3501 0.3506 0.3512 0.3517 0.3522 
1.031 0.3528 0.3533 0.3538 0.3544 0.3549 0.3554 0.3559 0.3565 0.3570 0.3575 
1.032 0.3580 0.3586 0.3591 0.3596 0.3601 0.3607 0.3612 0.3617 0.3622 0.3627 
1.033 0.3632 0.3637 0.3643 0.3648 0.3653 0.3658 0.3663 0.3668 0.3673 0.3678 
1.034 0.3683 0.3688 0.3693 0.3698 0.3703 0.3708 0.3713 0.3718 0.3723 0.3728 
1.035 0.3733 0.3738 0.3743 0.3748 0.3753 0.3758 0.3763 0.3768 0.3773 0.3778 
1.036 0.3782 0.3787 0.3792 0.3797 0.3802 0.3807 ~—-- 0.3812 0.3816 0.3821 0.3826 
1.037 0.3831 0.3836 0.3840 0.3845 0.3850 0.3855 0.3859 0.3864 0.3869 0.3874 
1.038 0.8878 0.3883 0.3888 0.3892 0.3897 0.3902 0.3906 0.3911 0.3916 0.3920 
1.039 0.3925 0.3930 0.3934 0.3939 0.3944 0.3948 0.3953 0.3957 0.3962 0.3967 
1.040 0.3971 0.3976 0.3980 0.3985 0.3989 0.3994 0.3998 0.4003 0.4007 0.4012 
1.041 0.4016 0.4021 0.4025 0.4030 0.4034 0.4039 0.4043 0.4048 0.4052 0.4057 
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0.01962 
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0.06820 
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TABLE 1 (Continued) 


OF FLOW OF AIR AND DIATOMIC 
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‘7151 
. 7300 
7441 
7573 
7699 
7817 
-7930 
8036 
8137 
8233 
8324 
8410 
8493 
8571 
8645 
8716 
8784 
8849 
8910 
8968 
9024 
9078 
9128 
9177 
9223 
9267 
9309 
9349 
9387 
9424 
9459 
). 9492 
9523 
9553 
9582 
“9609 
9635 
9660 
9683 
9706 
9727 
9747 
9766 
9784 
9801 
9817 
9833 
9847 
9861 
.9873 
9885 
9897 
). 9907 
9917 
9926 
9935 
9943 
9950 


5 


.4083 
.4127 
.4170 
.4212 
4254 
4295 
.4336 
.4376 
.4588 
.4939 
. 5254 
. 5539 
. 5800 
.6041 
6263 
.6469 
.6661 
. 6840 
. 7008 
. 7166 
7315 
. 7454 
. 7586 
7711 
. 7829 
. 7941 
8046 
.8147 
, 8242 
. 8333 
.8419 
.8501 
8579 
. 8653 
8723 
.8791 
. 8855 
8916 
.8974 
. 9030 
. 9083 
. 9133 
.9181 
. 9227 
.9271 
.9313 
. 9353 
. 9391 
.9427 
. 9462 
.9495 
. 9526 
. 9556 
. 9585 
. 9612 
. 9638 
. 9662 
. 9686 
. 9708 
.9729 
. 9749 
. 9768 
. 9786 
. 9803 
. 9819 
. 9834 
. 9848 
. 9862 
. 9875 
. 9887 
. 9898 
. 9908 
.9918 
. 9927 
. 9936 
. 9944 
.9951 


0.4087 
0.4131 
0.4174 
0.4216 
0.4258 
0.4299 
0.4340 
0.4380 
0.4625 
0.4972 
0.5283 
0.5566 
0. 5825 
0.6064 
0.6284 
0.6489 
0.6679 
0.6858 
0.7025 
0.7181 
0.7829 
0.7468 
0.7599 
0.7723 
0.7840 
0.7951 
0.8057 
0.8157 
0.8251 
0.8342 
0. 8427 
0.8509 
0.8586 
0. 8660 
0.8730 
0.8797 
0.8861 
0. 8922 
0.8980 
0.9035 
0.9088 
0.9138 
0.9186 
0.9232 
0.9276 
0.9317 
0.9357 
0.9395 


. 9665 

. 9688 

.9710 

.9731 
0.9751 
0.9770 
0.9788 
0.9805 
0.9821 
0.9836 
0.9850 
0.9863 
0.9876 
0.9888 
0.9899 
0.9909 
0.9919 
0.9928 
0.9937 
0.9944 
0.9952 


0 
0 
0 
0 
0.9640 
0 
0 
0 
0 


ooo sososs Sesss999c9099 


— 


oooco 


GASES 
7 8 
4092 0.4096 
4135 0.4140 
4178 0.4183 
4221 0.4225 
4262 0.4266 
4303 0.4307 
4344 0.4348 
4384 0.4388 
4661 0.4697 
5004 0. 5037 
5313 0.5342 
5593 0. 5620 
5850 0.5875 
6086 0.6109 
.6305 0.6326 
6508 0.6528 
6698 0.6716 
6875 0.6892 
7041 0.7057 
7197 0.7212 
7343 0.7357 
7481 0.7495 
7612 0.7624 
7735 0.7747 
7852 0.7863 
7962 0.7973 
8067 0.8077 
8166 0.8176 
8261 0.8270 
8350 0.8359 
8436 0.8444 
8517 0.8524 
8594 0.8601 
8667 0.8674 
8737 0.8744 
8804 0.8810 
8867 0.8873 
8928 0.8934 
8985 0.8991 
9041 0.9046 
9093 0.9098 
9143 0.9148 
9191 0.9196 
9236 0.9241 
9280 0.9284 
9321 0.9325 
9361 0.9365 
9399 0.9402 
9434 0.9438 
9469 0.9472 
9501 0.9505 
9533 0.9536 
9562 0.9565 
9590 0.9593 
9617 0.9620 
9643 0.9645 
9667 0.9670 
9690 0.9693 
9712 0.9714 
9733 0.9735 
9753 0.9755 
9772 0.9773 
9789 0.9791 
9806 0.9808 
9822 0.9824 
9837 0.9839 
9851 0.9853 
9865 0.9866 
9877 0.9878 
9889 0.9890 
9900 0.9901 
9910 0.9911 
9920 0.9921 
9929 0.9930 
9937 0.9938 
9945 0.9946 
9952 0.9953 
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9 
0.4101 
0.4144 
0.4187 
0.4229 
0.4271 
0.4312 
0.4352 
0.4392 
0.4733 
0.5069 
0.5371 
0. 5646 
O 5899 
0.6131 
0.6347 
0.6547 
0.6734 
0.6909 
0.7073 
0.7227 
0.7372 
0.7508 
7637 
7759 
7874 
7984 
8087 
8186 
8279 
8368 
8452 
8532 
8609 
8681 
8751 
8817 
. 8880 
. 8940 
. 8997 
. 9051 
. 9103 
.9153 
. 9200 
.9245 
. 9288 
. 9329 
). 9369 
. 9406 
0.9441 
0.9475 
0.9508 
0.9539 
0.9568 
0.9596 
0. 9622 
0.9648 
0.9672 
0.9695 
0.9716 
0.9737 
0.9757 
0.9775 
0.9793 
0.9809 
0.9825 
0.9840 
0.9854 
0.9867 
0.9880 
0.9891 
0.9902 
0.9912 
0.9922 
0 
0 
0 


S99 99999 S9D99 99999999999 


oo 


.9931 
. 9939 
. 9947 
0.9954 





(Continued on page 312) 
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TABLE 1 (Continued) 

Pu/pe tes tm 2 3 4 5 i Bn oe 8 Ga 
1.74 0.9954 0.9955 0.9956 0.9956 0.9957 0.9957 0.9958 0.9959 0.9959 0.9960 
1.75 0.9961 0.9961 0.9962 0.9962 0.9963 0.9964 0.9964 0.9965 0.9965 0.9966 
1.76 0.9966 0.9967 0.9968 0.9968 0.9969 0.9969 0.9970 0.9970 0.9971 0.9971 
1.77 0.9972 0.9972 0.9973 0.9973 0.9974 0.9974 0.9975 0.9975 0.9976 0.9976 
1.78 0.9977 0.9977 0.9978 0.9978 0.9978 0.9979 0.9979 0.9980 0.9980 0.9981 
1.79 0.9981 0.9981 0.9982 0.9982 0.9983 0.9983 0.9983 0.9984 0.9984 0.9984 
1.80 0.9985 0.9985 0.9985 0.9986 0.9986 0.9986 0.9987 0.9987 0.9987 0.9988 
1.81 0.9988 0.9988 0.9989 0.9989 0.9989 0.9990 0.9990 0.9990 0.9990 0.9991 
1.82 0.9991 0.9991 0.9992 0.9992 0.9992 0.9992 0.9993 0.9993 0.9993 0.9993 
1.83 0.9994 0.9994 0.9994 0.9994 0.9994 0.9995 0.9995 0.9995 0.9995 0.9995 
1.84 0.9996 0.9996 0.9996 0.9996 0.9996 0.9996 0.9997 0.9997 0.9997 0.9997 
1.85 0.9997 0.9997 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998 0.9998 
1.86 0.9999 0.9999 0.9999 0.9999 0.9999 0.9999 0.9999 0.9999 0.9999 0.9999 
1.87 0.9999 0.9999 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 





N = 1.000 for all higher pressure ratios. 


PROPOSED FORMULA 


In place of separate formulas for subcritical and 
critical flow, it is proposed to substitute the single 
general formula 


w= KCAp,,N/ VTi, 


which is merely the formula for critical flow multiplied 
by N, the tabulated restriction factor. The value of NV 
varies from zero to unity for subcritical flow and is unity 
for critical flow. 

This formula was introduced into the General Elec- 
tric Company some years ago by the late V. Petrovsky 





0 





4S te 


Fic. 1. Values of restriction factor N for different values of 
isentropic exponent k. 


and has been found to be convenient and timesaving for 
both air and steam calculations. Values of N may be 
read from curves such as are shown in Fig. 1. These 
are reasonably satisfactory but have the usual disad- 
vantages attendant upon purely graphical presentation. 
The formula is used so frequently and the advantages of 
a table over a chart are so manifest that an accurate 
table for air and diatomic gases has been prepared, 
based, for the sake of consistency, on fundamental con- 
stants which are essentially those adopted in the 
A.S.M.E. codes for flow measurement! and for the test- 
ing of air and gas compressors. These codes have been 
extensively used in the aeronautical industry, and it 
appears best to standardize on values already in current 


use. 
DERIVATION OF FORMULA 


The general derivation of the flow formula is well 
known, and the details will not be emphasized here. 
The case considered is that of the flow of air or of a 
diatomic gas through a nozzle or orifice of any kind 
with a velocity increase corresponding to the pressure 
drop. The pressure ratio is taken as the high-side fotal 
pressure divided by the low-side static pressuie, and this 
is obviously always greater than unity if there is any 
flow. The term “high-side” refers to entrance condi- 
tions and is somewhat more expressive than the words 
“inlet” or “‘initial’’ frequently used. The use of total 
pressure and total temperature at this section makes 
unnecessary the introduction of any approach-velocity 
correction. The usual assumption as to uniformity of 
flow over the cross section is made. 

The low-side static pressure is the pressure of the 
region into which the fluid is discharged. For sub- 
critical flow, this is the same as the pressure at area A 
(or throat), but for critical flow it may be, and usually 
is, less. Since it is only in the subcritical case that the 
low-side pressure affects the flow, it is not important in 
general (at least from the standpoint of flow measure- 
ment) to know the actual low-side pressure if it is 
known definitely that it is below the critical value. 

Considering first subcritical flow, the theoretical 
velocity head acquired in an expansion from inlet to 
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discharge pressure is equal to the change of enthalpy 
(otherwise called available energy, isentropic heat drop, 
or adiabatic heat drop), which is equal to the specific 
heat at constant pressure multiplied by the isentropic 
temperature drop. 


V3?/2g = Cp71.X/(X + 1) 





or 
Vs = V 2gcpTy,X/(X + 1) (1) 
The flow per unit of effective area is 
w/CA = V3/v; (2) 


Assuming, as usual, an isentropic process (pv* = 
constant), 








v3 = v1 7)/* = 14,(X a 1)1/@ = (3) 

so that 
wes eet BB acdat viaey 
CA a ni? (X _ 1)@ + D/@ —») ) 


It is apparent that w/CA is a maximum when the 
function 


o(X) = X/(X + 1)F + 0/4 -Y (5) 
is a maximum. Differentiating this expression with 


respect to X and equating the derivative to zero, the 
critical value of X is found to be 








X, = (k — 1)/2 (6) 
For the critical pressure ratio, then, 
We _ 267 X, a (7) 
CA ny? (X,+ 1)e+0/e-09- 


Dividing Eq. (4) by Eq. (7), the formula for the re- 
striction factor is obtained. 


N= poh = i SD lata 1) (8) 
we X, (Xx + 1) 


For a constant value of k, this expression is a function 
of X only. It may therefore be tabulated for different 


values of r. 
The flow for any pressure ratio may then be calcu- 


lated from the formula ° 


Bia y yet xX, (9) 
CA Pd (X, + 1)@ +)/— —D 








which is reduced to the customary simple form of the 


critical flow formula by using the relation 
Cf —c=R 
or 
co{1 — (1/k)] = prta/ Ts 


so that 
Me = (CoT11/Pid((R — 1)/R] (10) 


TABLE 2 
Values of K for Use in General Flow Formula 








Pu Tu A w K 
Ibs. per sec. 0.53033 


i 
sq.in. beg per min. 31.820 
°F. absolute libs. per hour 1909.2 
Ibs. per sec. 76.368 
sq.ft. {lbs. per min. 4582.1 
Pounds per sq.in. Ibs. per hour 274,920. 
absolute Ibs. persec. 0.39529" 
|sq.in. {Ibs. per min. 23.717 
~~. absolute Ibs. per hour 1423.0 
Ibs. persec. 56.921 
sate Ibs. per min. 3415.3 
Ibs. per hour 204,920. 
( (Ibs. persec. 0.26048 
|sq.in. (Ibs. per min. 15.629 
7. absolute} Ibs. per hour 937.74 
(lbs. persec. 37.510 


sq.ft. libs per min. 2250.6 
Ibs. per hour 135,030. 
(Ibs. persec. 0.19415 

sq.in. {lbs. per min. 11.649 
(Ibs. per hour 698.95 


Inches Hg abso- 
lute 


{ 
( 
‘. absolute 
| 
( 
if 





Ibs. persec. 27.958 
sq.ft. {lbs. per min. 1677.5 
Ibs. per hour 100,650. 
( ea persec. 0.0036829 
|sq.in. ;Ibs. per min. 0.22097 
oR. absolute] (Ibs. per hour 13.258 
Ibs. persec. 0.53033 
pane Ibs. per min. 31.820 
Pounds per sq.ft. Ibs. per hour 1909.2 
absolute ( Ibs. per sec. 0.0027450 
{sq.in. (lbs. per min. 0.16470 
°C. eae. Ibs. perhour 9.8821 
Ibs. persec. 0.39529 
lsq.ft. (Ibs. per min. 23.717 
| Ibs. per hour 1423.0 


Substituting Eqs. (6) and (10) in Eq. (9), the formula 
becomes 


w/CA = Kpi,.N/V Tu (11) 


K=[ k? (, 2 ae (12) 
k—1\k+1 c 


with 





For k = 1.39470, 
K = 1.2855 Vg/c, 


and 
w/CA = 1.2855 V(g/e)(buN/VTu) (13) 


which is the general formula for flow permitting the use 
of any consistent units. 

It is probably desirable to assume some particular 
system of units and reduce the formula still further in 
order to eliminate as far as possible the amount of calcu- 
lation required in routine computations. There is, 
however, no general agreement on standard units to be 
used for w, A, p, and T, and Table 2 has been prepared 
to give the value of K corresponding to all the different 
combinations of units likely to be used in practice. The 
working formula for flow for any given values of C, A 
Pi» Ti, and r then becomes 
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w= KCAp,,N/ VTi, 
with N taken from Table 1 and K from Table 2. 


SELECTION OF Basic CONSTANTS 


The basic physical constants used in preparing the 
tables are essentially those adopted in A.S.M.E. test 
codes for measurement of flow and for compressor test- 
ing. In particular, the following constant values have 
been used directly or are implied. 


po = 14.70 Ibs. per sq.in. absolute 

T) = 527.6°F. absolute (68°F.) 

yo = 0.075 Ib. per cu.ft. (exact value by definition) 
(k — 1)/k = 0.283 (exact value by definition) 

k = 1.39470 

C. 0.242967 B.t.u. per Ib. per °F. 

g = 32.1740 ft. per sec. per sec. 


The “‘long”’ values of k and c, were used in computing 
the constant A in order to ensure exact consistency 
with the selected value of (k — 1)/k, although even 
here, where the computation needs to be made only once 
and for all, this is a rather academic proceeding. Nor- 
mally, these two values would be rounded off to three or 
four significant figures. 

The reasons for selection of these constants are ex- 
plained in the A.S.M.E. codes and in the references 
there given. Since it may not be convenient in all 
cases to consult these references, the reasons will be 
briefly summarized here. 

First, ‘‘Normal Air”’ is defined as air with an average 
amount of CO, and sufficient moisture so that its specific 
weight at the selected standard conditions—p) = 14.70 
Ibs. per sq.in. absolute and Ty) = 527.6°F. absolute 
(68°F.)—is exactly yo = 0.075 lb. per cu.ft. These 
values were selected chiefly because they (or values sub- 
stantially equal to them) were already being used as 
standards in allied fields, especially in heating and 
ventilating engineering. 

By combination of the two fundamental thermo- 
dynamic relations 


Cp — & = c,[1 — (1/k)] = R (14) 
and 
pv = RT (15) 
the relation connecting c,, k, and v (or y) is obtained 
ell —(1/k)] = po/T=p/yT (16) 


Investigation of published data on the specific heat of 
air at temperatures of usual interest (0-200°F.), and 
containing about the same amount of moisture as 
Normal Air, led to the selection of 0.243 B.t.u. per Ib. 
per °F. as a reasonable average value. Substituting 
this in Eq. (16) and using on the right-hand side the 
values applicable to Normal Air, it was found that 


(k — 1)/k = 1 — (1/k) = 0.28296 





TICAL SCIENCES—JUNE, 1946 


Since (k — 1)/k is a function more commonly used in 
thermodynamic calculations than & itself and since the 
value derived from Eq. (16) can be only an average 
value at best, this figure was rounded off to 0.283, which 
for the purpose of establishing a consistent set of values 
for thermodynamic calculations was taken as exact. 
The corresponding, value of c, then became 0.242967, 
which would be rounded off to 0.243 for usual calcula- 
tions but must be retained in full when exact consistency 
is required. Similarly, the corresponding consistent 
value of the ratio of specific heats isk = 1.39470. 

The ratio (k — 1)/k is selected as a basic value be- 
cause many important thermodynamic formulas, in- 
cluding formulas for the flow of compressible fluids,’ 
may be expressed in terms of a certain function 


(br/ps)* -9/* — 1 


which is sometimes known as ‘‘adiabatic factor’’ but is 
probably better designated as ‘“‘isentropic’’ factor. 
Tables of this function are now available in various pub- 
lications, as well as in A.S.M.E. codes. It has been ex- 
tensively used in the aeronautical industry, usually be- 
ing denoted by the symbol, Y, which was originally 
given to it.° With increasing use, however, the 
American Standards Association has decreed that, in 
order to avoid certain conflicts, it shali be denoted by 
the symbol X rather than Y, and it is so designated in 
most recent publications. The symbol X is therefore 
adopted in this paper. 


UsE OF FORMULA 


The formula is given in terms of ;, and 7),, the total 
pressure and temperature, respectively, on the upstream 
side of the restriction. In most cases the reading of the 
usual thermometer or thermocouple is nearly enough 
equal to the total temperature, but for extremely high 
velocities the reading is usually intermediate between 
the true static and the true total temperature and a 
correction must be applied. ‘This is true, of course, no 
matter what formula is used. 

The total pressure is the pressure that would be 
shown by an impact tube pointing directly upstream. 
Sometimes a static rather than a total pressure is 
measured. The static pressute and temperature may 
be used in the formula if the usual correction for velocity 
of approachis made. This correction is often negligibly 
small. 

In making flow measurements it is usually considered 
better practice to measure the differential 


Ap; = pu — pr 


directly, together with either one or both of the gage 
pressures, rather than to depend only on the latter. 
The reduction to common units and calculation of the 
pressure ratio is combined into one operation by the 
formula 
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CALCULATION OF FLOW OF 


Pu _ 4 K’' Ap, B 1 

po p2 1 — (K’Ap,/pu) 
Values of K’ are given in Table 3 for the cases most 
likely tooccur. The constants there given assume that, 
in accordance with usual practice, the readings of mer- 
cury columns have been reduced to equivalent heights 
at 32°F. and that the readings of water columns have 
been reduced to, or are considered closely enough equal 
to, equivalent heights at 60°F. 

The formula gives what is called “‘gravimetric’’ flow. 
This term indicates merely a flow that the practical 
engineer is accustomed to express in pounds per second, 
pounds per minute, or the like and takes no stand on the 
controversial issue as to whether it should be denoted 
“mass” flow or “weight” flow. The corresponding 
volumetric flow at any specified pressure and tempera- 
ture is obtained by multiplying by the specific volume 
at that pressure and temperature. 

In certain problems the formula is used in the form 


N= (Vv T,,w/A)/KCpi 





as, for example, when checking the capacity of a restric- 
tion or duct for a specified flow. If the calculated value 
of N is greater than unity, it is impossible to pass the 
specified flow with the specified pressure and tempera- 
ture. With N exactly equal to unity, the velocity is 
sonic. A value of N less than unity indicates the frac- 
tion of the total possible capacity being utilized. 


CALCULATION OF TABLES 


Most of the values of NV have been calculated using 
seven-place logarithmic tables and a _ calculating 
machine. For the Jow-pressure ratios, where the num- 
ber of significant figures given in seven-place logarithmic 
tables is rather small, the 16-place table of natural 
logarithms prepared by the Federal Works Agency 
under the sponsorship of the National Bureau of 
Standards* was used. In some cases the logarithms 
were themselves calculated directly from the appropri- 
ate infinite series. 

Since the function NV changes rapidly at low-pressure 
ratios, the increments of pressure ratio have been chosen 
with a view to making the differences between the 
tabulated values small enough for convenient interpola- 
tion throughout the entire table. However, interpola- 
tion should not be necessary in most practical cases. 
The table was originally prepared with the values of NV 
carried out to the nearest fifth significant figure, which 
for the greater portion of the table is equivalent to five 
decimal places. However, it seemed that this number 
of places was greater than was needed in the great 
majority of practical cases and, in fact, would seldom be 
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TABLE 3 
Values of K’ for Use in Calculating Pressure Ratio 


Pu/p2 = 1 + (K’Ap,/p2) = (Ran 








Pur P2 Ap, e 
(Inches H,Of 0.07347 
Inches Hg absolute* {Inches Hg* 1.000 
trenaas per sq.in. 2.036 
{Inches H,O+ 0.03609 
Pounds per sq.in. absolute; Inches Hg* 0.4912 
Pounds per sq.in. 1.000 
* At 32°F. 
+ At 60°F. 


justifiable when the accuracy of the usual test and other 
data available are taken into consideration. The values 
here given have therefore been rounded off to the 
nearest fourth significant figure, with the incidental 
advantage of facilitating interpolation if this is re- 
quired. 

It is shown in reference 2 (formula D) that for low- 
pressure ratios the flow is closely proportional to 


prizcas 
V poAp, 
Equating this closely approximate formula to the 


formula here used and solving for JN, it is found that for 
low-pressure ratios N is proportional to 


V bs Ap./ pis 


If all pressures are expressed in the same units, the 
constant of proportionality varies from 2.06806 at a 
pressure ratio of unity to 2.06416 at a ratio of 1.0500, 
with intermediate values nearly in direct proportion. 
Occasionally, this relation may be of use in calculating 
values of NV not tabulated directly. Since N is a func- 
tion of pressure ratio only, any convenient values of 
absolute pressure may be assumed in making the calcu- 
lation. 
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Coupled Surfaces 


J. B. KENDRICK* 
California Institute of Technology 


INTRODUCTION 


HE OBJECT OF THIS INVESTIGATION is to describe 
7. new type of airplane control surface that is many 
times more powerful, as a stabilizing or damping means, 
than a fixed surface of equal area. This increase in 
effectiveness is accomplished in the following manner: 

As shown in Fig. 1, the lateral force F (1) produced by 
angle of yaw y (2) causes the vertical surface to tip (3) 
about its base. Coupling link (4) between the rudder 
and a fixed point (preferably inside the fuselage) causes 
the rudder to deflect through angle 8 (5). The rudder 
deflection acts in the proper direction to reduce the 
angle of yaw y. This procedure has converted the con- 
trol surface into an automatic servo-control. 











General view of coupled fin and rudder 


Fic. 1. 


Although Fig. 1 illustrates a vertical surface, this 
principle is not limited in usefulness to vertical tails 
alone. Any control surface that is hinged at the trail- 
ing edge of an airfoil may be adapted to an automatic 
servo-control in two steps, as follows: 

(1) Hinge the forward portion of the airfoil upon a 
fore-and-aft axis. This provides the “‘error-smeller’’ 
element of the servo-control, because the forward part 
will tip whenever the surface is deflected from a trimmed 
position. 

(2) Couple the rearward or flap portion to a fixed 
point on the airplane in such a manner that movement 
of the forward portion on its hinge will cause the flap 
to move in a direction to provide aerodynamic correc- 
tion for the error in trim. In this way, the coupling 
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linkage serves to actuate the flap and also to stabilize 
the forward part of the surface. 

When the foregoing procedure is carried out accord- 
ing to certain rules that are to be given, the resulting 
control surface is sensitive to changes in angle of attack 
and acts automatically to reduce the error with a force 
proportional to the error.’ (In this respect it functions 
in the same manner as an angle-of-attack-type auto- 
pilct.) As a free control, its effect upon the airplane 
is the same as a conventional surface of much larger area. 
It may be adapted to manual control in the normal 
manner, but one should note that the maximum lift 
or moment will be no greater than a conventional sur- 
face, even though the slope of the lift curve is many 
times as great. 

Since any servomechanism is susceptible to a tendency 
to oscillate unless properly damped, this subject has 
been carefully studied, and stability criteria are pro- 
vided for use in design. The resulting surface will re- 
turn promptly to its original setting when deflected from 
a trimmed position in an airstream. 

The static and dynamic stability will now be developed 
for the case of a vertical tail surface. This develop- 
ment will likewise serve for horizontal tails and, with 
slight modifications, may be adapted to the case of 
ailerons. 


DEFINITION OF TERMS 


total mass of tail 

m, = mass of rudder 

c.g. A to rudder hinge 

x = cg.A toc.g. complete tail 
rudder hinge to rudder c.g. 
c.g. A to tip hinge 

tip hinge to c.g. complete tail 
tip hinge to c.g. rudder 


J 
= 
= 

~ 


8 
° 
I 


83 
tou tl 


2, 


t Another arrangement of all-movable surfaces is described in 
reference 3, where the forward surface is hinged on an axis 
parallel to the rearward flap surface. Although this method of 
coupling provides a considerable increase in dC,/da, the ar- 
rangement of perpendicular hinges described herein is even more 
powerful, as shown by the following comparison of the effective- 
ness factor [1 + r(8/W)] for the two methods: 





Static Dynamic Source 
Method of Coupling Condition Condition of Data 
: theory >2 chai erat 
Parallel hinges he ie. 1.15! Reference 3 
fa - theory >8 >6 Figs. 3 and 4 
Perpendicular hinges { test 76 2.02 Table 2 





1 Highest value reported. Larger values should be feasible. 

? The dynamic stability theory was not completed when these 
tests were made. See “Correlation of Dynamic Stability Theory 
with Test.” 
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: . zi tip z 
R = opting © 0<R<1;G=-R 
ZI x 
S = total tail area 
S, = rudder area 
C = total tail chord 
C, = average rudder chord 
% = cg. A to '/< complete tail 
t, = 1/,c complete tail toc.p. rudder force due to 8B 
€ = tip hinge toc.p. complete tail 
C dC, a dC, dCi da’ ‘*) 
sag da ’ “ dp ’ ar da ’ aT dp 
z, 
= 1 = 
Ti T ( a : ) 
CraSrCr CpS-C; Cimgg C 
he = a hg = Oe 1 + ue) 
CraSé CraS2 Creat 2 
MyXr ' X—-x+x,* 7p IC, (> 
M = —; (1+f*) = ~———; Cag = —— 
mx x dp 


a and } = radii of gyration of rudder (about a vertical axis) 
and the complete tail (about a horizontal axis), respec- 
tively, through their own c.g.’s 


g ; +, Se 
H = 3 he — rR) + [om + = - 26) ] + 
§ x z 


A 
a 


G— MA +f*) |] 1 — A(R — ha) 
~~" Se 2 
J = G+ uf ant 26%) - E — M(1 +7] 
K = —[n(R — ha) + (hg — 7R)]; L = ——— 
3 Mx 
y = angle of yaw 
8 = rudder deflection 
a = angle of attack of airfoil 
q = dynamic pressure = (p/2)U? 
p = air density 
U = air velocity 


StaTiIc STABILITY OF COUPLED FIN AND RUDDER 


The resistance to sideslip dC,/dy is equal to the sum 
of the lift-due-to-angle-of-attack dC,,/da and lift-due-to- 
rudder-deflection (dC,/da)[(da’/d8)(dB/dy)|, where 
da’/dB is the change in effective angle of attack of the 
tail per degree of rudder deflection, normally denoted 
by 7; and d8/dy is the rudder deflection per degree of 
yaw, which may be written 8/y for small deflections. 
The sum of these effects is 


dC, oz | s) : 
en wt cae 2 = 1) 
ie eae 


The dC,/da of the normal tail is thus augmented by the 
factor 7(8/~)(dC,/da) produced by the forced rudder 
deflection. The coupled fin and rudder are therefore 
[1 + 7(B/yp)] time as effective as a fixed vertical surface 
of the same area. The value of 8/y may now be derived 
in terms of the geometry. (The symbols are given 
above.) 

As shown in Fig. 2, aerodynamic tipping moment = 
link restraining moment 


(dC,/da)(y + 78)gSz = —F(y,/l)2, (2) 


and rudder hinge moment = link restraining moment 


(Chav + CrgB)gS,C, = — F(y,/l)x;, (2a) 
Dividing Eq. (2a) by Eq. (2) and solving for 8/y, 
B/pP = (R — hy)/(hg — 7R) (3) 


As the coupling ratio R increases, the effectiveness 
{1 + 7(8/yp)] likewise increases, until at some value 
R, corresponding to hg — +R,, the denominator of Eq. 
(3) is equal to zero and [1 + 7(8/y)], approaches in- 
finity. At larger values than R, the coupled surface is 
unstable. R, and its effect on stability are determined 
as follows: 


i tuefie 1 + AR = Me). hg — rR, = 0 
y hg ~ TR 
therefore 
R, = hg T 

R ha R.,e 
Mapa ay see ey 

eo ae Se. © de zien Op 
v - 1— (R/R,) (R;/R) -— 1 


(3a) 


Optimum values of the effectiveness [1 + 7(8/y)] 
may be obtained in the following manner: 
(A) Make the numerator as large as possible, and 


positive. 
(1) Large 7 indicates that a wide chord rudder 
should be used. 


(2) Large R may be obtained with wide chord 
rudder on a low aspect ratio tail. How- 
ever, R is really determined by B(2) be- 
low. 

(3) h, should be as small as consistent with 
large hg under B(1). 

(B) Make the denominator as small as feasible, and 
positive. 

(1) Large tg = CigC,S,/C,.2S indicates: 

(a) Low aspect ratio, to give small 2 and, 
incidentally, low C,,. 

(b) Large rudder chord C, and area 
ratio S,/S. 

(c) High hinge moment coefficient 
slope, Cyg, suggests a rudder with 
small aerodynamic balance, and 
with cusp trailing edge. Surface 
contours should be carefully con- 
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Fic. 2. Coupling link geometry. 
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trolled to avoid nonlinearity of 3 Q 
Cg near neutral. I 14th n+ PEt 0 
(2) 7rR should be adjusted to slightly less e-! = 
thanhg. Aerodynamic and manufac- ‘ | Theory 
turing variations will prove the de- Test x 
ciding factor in the minimum value to é 
be used for (kg — +R), as shown in | 
) 
Eq. (8a). . 
Allowing a margin of say 10 per cent from R, to cover e | f 
aerodynamic and manufacturing variations, one ob- uj. 
tains the following value of effectiveness for normal sur- 7? . 
face proportions, assuming r = 0.65 and h,/hg = 1/2. | | | d 
1 + 7(8/p¥) = 1 + {7[(0.9hg/7) — hq) /0.10hg} = ai | | dt 
1 + [0.9 — (0.65/2)]0.10 = 6.8 : NK | | | 
Ler : ¥ xy S ] 
A decided improvement in static stability may therefore G a d 
be effected with surfaces coupled in this manner. DN see seen dl 
Excellent agreement is obtained between the static 7 
stability theory and test results, as shown in Fig. 3. Ratio ¥, 
The experimental points were obtained from free con- 40 2 | el, - a d 
trol wind-tunnel tests of the coupled tail at various reer eae, A eae i 
angles of yaw, as shown by the te Fig. 4. 1G.3. Static stability characteristics of coupled fin and rudder, ( 
Rudder deflection tests (as shown dotted in Fig. 4), d 
when cross-plotted for various linkage ratios, gave the a £ vee row. o Coupling Lin a dt 
same results as the coupled tail. > ae 1 o wt at, 
E ag aeeeers 
DyNAMIC STABILITY OF COUPLED SURFACES ‘ sae ae a eo 
ja & AF Ge Ao 
Equations of Motion is x Ke ane Te ae 
The motion of the surfaces may be expressed in terms * 8 ae, V, ‘ / ae : |, i ai 
of the angle of yaw y and the rudder deflection 6. It @ 4 Ke UW) 7 Va fyi dP 
is now desired to express ¥ and £ in terms of equations ae ; Hy - VY vO ne d 
of motion with respect to time, in order to indicate the ff 2 A Lx A~ 
dynamic behavior of the coupled surfaces. Further- y? ge /, Vf a 
more, it is desired to develop a stability criterion that Ja’. Y ghd 
shows the relative importance of the various aerody- 44 5g sin 
namic and mechanical factors, and their proper rela- ai po Angld of Yaw\ ¥ - Degrees str 
tionship to insure satisfactory stability. This may oO. +P 7 4 8 i ‘ 2 
be done with the aid of Lagrange’s equations of mo- Al 
tion.! 
(d/dt)(OT'/0gx) — (OT / qx) =. Qx (4) A} 
In Fig. 5, any small element P; of mass m; and Carte- 
sian coordinates x;, y;, 2; may be expressed in terms of 
y and 8, the generalized coordinates. 
x, and z, will be constant for small values of y and 8. 
yi = (xo + x,)¥ — R28 (for fin) ‘ 
Vi = (Xo + xv — (Re, — x;)6 (for rudder) 
9: = (xo + x)~ — Rz,B (for fin) 
9, = (% + x)) — (Re, - x)B (for rudder) 
%,=2,=0 
my 





Fic. 5.. Geometry of coupled fin and rudder. 


The total kinetic energy 7 and its derivatives are used 
to obtain the inertia moments of the surfaces. a , ae 
T = '/2>>m,[(xo + xb — RziB]? + 


K.E. = T = 1/,> 0m (%2 + 92 + 2) + 1/,>0m wv? 1/5°m 4|(x0 +- x.) ~ (Rs, — x ;)6)? Det 
i=1 


i=! 











a 
_ad 


_ rudder, 





x ;)B]? 
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4 = 20m [2(xo + x1)? ¥ — 2x0 + x1)ReiB] + 
(fin) 
"/ 2m ;[2(xo + x;)*y - 2(x0 + x;)(Rz; 


(rudder) 


wi x:)B| 


i = '/20m[—2(x0 + xR + 2R*228) + 
Cc ’ 
1/2>9m ;[—2(xo + x1)(Rz; = xv + 
2(Rz; — x;)°6] 


a( 2) Yim [(xo + x:)*b — (xo + x) R26] + 
dt\ Op i 
Lim l (xo + x1)*P — (wo + x1)(Rei — x,)B| 


df oT : , 
427) = Dm,[—(xo + x) Regd + R278] + 
dt\ OB i 

dom [— (x0 + x:)(Rei — xp + (Rei — x :)?B] 


d af) . z os 
—} = mxy — mxzRB + m,x,(xo + x,*)B 
Abe 4 2 


(2t) _ 
dt\ op 
m,X,(xo + x,*)b — 2m,x,2,RB + m,x,2*B 
Note that 07 /d0y = O7/06 = 0. 
The potential energy terms Q, are represented by the 


aerodynamic moments about the airplane c.g., and the 
rudder hinge: 


—m x2zRp + m2**R°B + (4a) 


1P.E. 2 - Xi Ii Zi 
“~=0,= B (x, S54 vy, 2% 4 2, 281) - 
dx i=l O”gx Ode Og: 
XY, Oy 
i Od: 
since X; = Z; = 0. Considering (1) as chordwise 


strips of the airfoil of width dz 
AY, = C.. | «w + 7B) + (bv + 7B) sf qgCdz_ (both) 
ay, = | (Cat + Crg8) + (Crab + Cu) 2 Jaca 

‘ (rudder) 





Oy,/ OW = (x) + x;) (both) 
Oy ,/0B = —Rz; (fin) 
—(Rz; — x;) (rudder) 


(dy;/0B = 


Moments about the airplane c.g. and moments about 
the rudder hinge are, respectively, 


Q = ECia0C| + 7B) + (W + 78) “x 
(xo + x,)dz = Cat SH (¥ + 78) + 
: 7 
(¥ + 7B ‘] 
v ) 7 
Qs = LCi | w + 78) + (Wt |x 


r 


qC(—Rz,)dz + E] Crh + Cag8) + (4b) 


(Chav + C8) ] qC,*dz 
= ~Cal + 78) + Wt ba |x 


gSzR + | (Crav + Cng8) + (Crab + Cre B) =| x 
gSC, 


Treatment of incremental terms, the values of which 
are often quite important to the final answer, will be: 


(1) Terms marked (_ )* are evolved from }>m,x? 
‘ 


and }°m,z?, the moments of inertia of the rudder about 
i 


its hinge and of the complete tail about the tip axis, 
which become m,x,?[1 + (a@?/x,?)] and m,z?[1 + (6?/s*)|, 
respectively. For brevity these terms will be marked 
( )* for identification. 

(2) The centroid of 8-produced aerodynamic forces 
is aft of the quarter-chord point, at # + #,.* Let 7 
designate the product r[{l + (#/®}) = r{l + 
Cmog(C/®)(1/Crar)) - 

(3) The ratio (x» + x,*)/x will be expressed as 

(xo — x + x,*)/x. 


(1 + f*), where f* = 


Substituting values from Eqs. (4a) and (4b) into Eq. (4), the equations of motion then become 


mxy — Caai( # + v) + mx*[M(1 + f*) — GIB — CuagSin( 2 B+ )= 0 


T 


t Z = * - 
mx*[M(1 + f*) — GW — CragS2(ha — R49 - v) 4 mae Om —2MG*%+M~ | on 
1 Zz x 


CragSt(hg — R)( #8 + z)= Se 3 


Denoting d/dt by the operator \ and the coefficients by a, b, c, and d, Eqs. (4c) have the form 
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(aad? + ard + ao)y = —(ded® + DA + bo) 8 (4d) 
(ccd? + ad + coy = —(d2d*? + dr + do)B (4e) 
Dividing Eq. (4d) by Eq. (4e) gives the following quartic equation in ): 


(@edz — bat2)d* + (dedi + aibe — boc, — dice)d* + (dado + aod, — deco — boc, + aid, — dics)? + 
(aydo + aod, — dio — bots1)A + (odo — boo) = 0 (5) 


or 
An + B+ O17? + DA +E=0 (5a) 
Solving for the coefficients A, B, C, D, and E | 























a = m x? be = mx?[M(1 + f*) — G] Q= mx?[M(1 + f*) eas 1 ame om + u(*“-2G*) | 
x Z | 
ae en a on ners . 
rie CragSt hs CrgSi ‘ bia Cra SE (y 48 rR): i. Coad SE (y A rR) = 
U/é U/z U/z z U/z z 
eo enese eeeapeen co = —CreqgS#(h, — R)* |do=—CragSE(hB — rR) = 
A = (m x*)?{G*® + M[(x,*/x) — 2G(z,/2)] — [G — MQ + f*)]?} = (mex*)*{ J} 
_ 5 2 5 * - 
p. -—Saee 2 (he — rR) + G*2+4+M (= ~ *) 4 E ~ M( +79] ee (R—- he) = 
U/é z x z £ 
— Cr .q9SE 2 
—— H 
He Oe 
C = B(U/%) + E(#/U)? = B(U/) (6) 
D = 2E(%/U) 
E = (CiagS#)*{ (2/2) [1(R — ha) + (hg — rR))} 
E = L*(m,x*)*{ (2/8)(hg — rR)[1 + 1(R — ha)/(hp — 7R))} 
E = Lima)? {= a rR)| (1 + 72) Coop © Al = Lm x)?K 
x y Ce. x v 
2 (hg — rR) + G4? + u(% - 26%) + E —M(1+ | E «~~ he) 
a — CraqgSé {EF whe Ax AL OR = a 
A mx?(U/2) G*? + M[(x,*/x) — 2G(z,/z)] — [G - MA + *)? ~ U/e J 
(7) 
C/A = L(H/J) + [L/(U/2))}(K/J) = L(H/J) (8) 
D/C = [2/(U/#)(E/C) = [2L/(U/2))(K/H) (9) 
E/C = [(E/A)/(C/A)] = L(K/H) (10) 
bios <a [\* + (B/A)d + (C/A)] [M+ (D/C) + (E/C)] = 0 
The quartic Eq. (5a) represents stable motion, pro- (5b) 
viding all the coefficients A, B, C, D, and E are positive iat ‘i 2s ; : 
aig mage presents the origina] equation accurately, pro- 
and Routh’s discriminant (cf. reference 1) BCD — viding Sis ghia. . where 





D? — B*E > 0. The latter condition may be reduced 
to the form 


B\? E/z\? H K 
—} ~—-—)>@er-—4(—}>0 (11) , 
(3) 42) se (*) () and 


Eq. (5a) may be solved approximately by factoring into 
two quadratics. 


_ md Hj pet é Lu 
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But according to Routh’s discriminant, H//J — 4(K/H) 
> 0 (Eq. (11)), which means that |Ai2| > 2/As, 4]. 
Therefore |di,2| >>|%s,4| and the factors of Eq. (5b) 
are satisfactory for the purpose of this investigation. 
Both of the above modes represent oscillatory motion, 
as indicated by the fact that C/A > (B/2A)? and 
E/C > (D/2C)?. 

The frequency and damping characteristics of the 
coupled tail may now be compared to those of a fixed 
tail of the same area, as shown in Table 1. It will be 
noted that both modes of the coupled tail may have 
greater damping and high frequency than the fixed tail. 
The degree of improvement depends upon the proper 
adjustment of the various terms comprising H, J, and 
K. 

Optimum Proportions for Dynamic Stability 


The damping criterion is represented by Routh’s 
discriminant, which in -this case reduces to H/J — 
4(K/H) (Eq. (11)). The ratios H/J and K/H are 
found to be the same as those that determine the 
damping and frequency of the lst and 2nd modes, re- 
spectively, in Table 1. In order to attain the maxi- 
mum damping and frequency of the 2nd or long mode, 
H/J — 4(K/H) is made small but positive, by keeping 
the value of 7 small, J and K large. 

The term H/ is made small in the following manner: 

(1) Make hg — rR small but safely +, with allow- 
ance for variations. 

(2) Make [G — M(1 + f*)][1 + (CGnog/Crar)(C/%) | 
small but safely +, with allowance for variations. 

(3) Make x,*/x < 2G(z,/z), by adjusting 2, and x,*. 

The term K is made large by the same rules as for 
static stability, see “Static Stability of Coupled Fin and 


M [(x,*/x) — 2G(z,/z)]. In this regard a compromise 
must be effected between this term in J and the same 
term in H, in order to obtain the desired small value of 
the ratio H/J.* 


Summary of Damping Criteria 


In the preceding section it was shown that certain 
terms within Eq. (6) are also of particular importance, 
because if these terms have the wrong sign the adverse 
effect on stability may be greater than the nominal 
value of the term. The overall stability may, therefore, 
be changed by any one of the following factors, which 
should be adjusted as specified, for proper stability: 


(1) A, B,C, Dand E> 0 
(2) (H/J) — 4(K/H) > 0 
(3) (hg — tR) > 0 
(4) (R—h,) > 0 


@) [o-wa+ypy](1+ S84) > 0 


Lat x. 


Correlation of Dynamic Stability Theory with Test 


Before the stability theory was developed in its final 
form, dynamic wind-tunnel tests were conducted on two 
different model tail surfaces, with many different bal- 
ance conditions. When the stability criteria were used 
to predict, by analysis, the behavior of these models in 
the dynamic tests, the results are as shown in Table 2. 
In eight out of nine cases, the above criteria gave the 
correct answer on stability. 


* Typical relative values of the above factors for a vertical tail, 
might be H/J — 4(K/H) = 15/1 — 4(50/15) = 1.7. In this 
case, the time to damp to '/, amplitude for the 2nd mode would 
be 1/(2 X 3.3) = '/sth as great as a fixed surface of equal area, 
and the natural frequency would be 1.8 times as great. The Ist 























Rudder.’” The term J is made large by adjusting the mode would have a much higher damping and natural frequency 
rudder mass balance and mass_ distribution, than the 2nd mode. 
_— “ene TABLE | Qitreiey HIME) Testes SHERCOT 0 
Coupled Tail vs. Fixed Tail 
€ry ——_—_—_————Coupled Tail———_—_—____—~ 
Ist neg 2nd Mode 
ie Tail short period long period 
i . B C D E 
Equation of motion v+e7 gi + Ly =0 M+ A+ = r ak 
L rc l A ey C 
Damping ratio é a ae VL Lo! vL V2 PNE- C vii 
2U/EVSL  2U/z 2ANC 2U/kENT 2CNE U/t NH 
Undamped natural frequency w» V/L V/ C/A = VLVH/I V/E/C = VLVK, H 
dil ree | VHT 2\/K/H 
£ fixed 
iia So eeeet VA/T VK/A 
SS ee a 
Ratio = & Ist mode vz — ie 
£ 2nd mode 4K 
“ lst mode H H 
R wo on|7 — 2-2 
_ wo 2nd mode J 4K 
Ratls tiene to d to 1/ litud Ist mode _ 0.693/(B/2A) “i 3(; K 1) «! 
€ to damp to 1/2 ampituc’ Ind mode  0.693/(D/2C)  2\H H 
—_— - 2nd mode L Rant x 1H 
1 TE aig gaat, eee, See ged 
Ratio time to damp to !/2 amplitude a ro Hl UeH~ 2K 
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TABLE 2 ae it 
Dynamic Stability Correlation ge 
Model ] — | = a. a 
peer senda | | | Prop. 5 
Cond. | S&S 2 Si | es Se Spoogug af) ep ves tg ( 
mz, Ibs. | 1.98 1.95 1.54, | 1.54 | 1.54 1.54 | 1.62 | 2.68 | 2.68 2.9 of 
m,, lbs. | 0.38 0.38 0.38 0.38 0.38 0.38 0.48 | 1.15 1.15 1.5 he 
x, in. 36 36 36 36 36 36 36 | 84 84 84 va 
Xr, in. 1.38 1.38 1.38 1.38 1.38 1.38 1.00 3.2 3.2 iF: 
z, in. 7.75 6.06 3.75 3.75 3.75 3.75 3.75 8.4 8.4 8.4 Ce 
zp, in. 4.13 413 | 4.13 | 4.13 | 4.13 | 4.13 | 400] 8.9 | 8.9 8.4 pr 
R, 0.115 0.0955 0.0722 0.0955 0.115 0.161 0.0722 0.167 0.143 0.235 an 
S, sq.in. 84 84 84 84 84 84 84 2.64sq.ft| 2.64 2.64 th 
S;, sq.in. 40 40 40 40 40 40 40 1.39 1.39 1.39 
C, in. 8.4 8.4 8.4 8.4 8.4 8.4 8.4 &e iuy seg 
z, in. 34 34 34 34 34 34 34 s2 | 82 82 
&, in. 4.0 4.0 4.0 4.0 4.0 4.0 4.0 6.0 6.0 |}: 6.0 
Cha — 0.004 |— 0.004 |— 0.004 |— 0.004 |— 0.004 |— 0.004 |— 0.004 |— 0.004 — 0.004 |— 0.004 
Crp — 0.010 |— 0.010 |— 0.010 |— 0.10 |— 0.010 j;— 0.010 |— 0.010 |— 0.008 — 0.008 |— 0.008 
CLa — 0.04 — 0.04 |— 0.04 |— 0.04 |— 0.04 |— 0.04 |— 0.04 |— 0.04 — 0.04 |— 0.04 (si 
T 0.65 0.65 0.65 0.65 0.65 0.65 0.65 0.67 0.67 0 67 pe 
he 0.047 0.047 0.047 0.047 0.047 0.047 0.047 0.088 0.088 0.088 y 
he 0.118 0.118 0.118 0.118 0.118 0.118 0.118 0.175 0.175 0.175 of 
M 0.00735} 0.0075} 0.0095) 0.0095} 0.0095) 0.0095) 0.0082 0.0164 0.0164 0.0164 
2/x 0.215 0.169 0.104 0.104 0.104 0.104 0.104 0.100 0.100 0.100 
Xe/% 0.0383 0.0383 0.0383 0.0383 0.0383 0.0383 0.0278} 0.038 0.038 0. 047 
2/2 0.533 0.68 1.10 1.10 1.10 1.10 1.07 1.06 1.06 1.80 
2/% 0.118 0.118 0.118 0.118 0.118 0.118 0.118 0.073 0.073 0.073 ine 
1 + (a?/x,?) 2.33 2.33 2.33 2.33 2.33 2.33 3.6 2.33 2.33 2.33 
1 + (b?/s?) 1.25 1.5 2.0 2.0 2.0 2.0 2.0 1.5 1.5 1.5 
err ae EN i. 
1 + (&,/2) 1.06 1.06 1.06 1.06 1.06 1.06 1.06 1.07 1.07 1.07 ( 
1 + f* 1.09 1.09 1.09 1.09 1.09 1.09 1.10 1.09 1.09 1.09 str 
(1) (H/J) — 4(K/H)| 16.4 12.4 een 2.6 2.1 —45 6.4 — 3.4 — 4.5 6.6 sit 
3 S S S U S U U S 
(2)he — rR 0.043 0.056} 0.071} 0.056] 0.043] 0.013} 0.071 | 0.063 0.079 | 0.018 of 
S S S S S ; S Ss Ss Ss oat 
(3) R — ha } 0.068 0.048 0.025 0.048 0.068 0.114 0.025 0.079 0.055 0.147 qu 
| ; S S Ss Ss gory Ss Ss S pe 
(4) G-— M(1+ f*) | 0.0168 0.0080} — 0.0029)}— 0.0004 0.0024 0.0063|/— 0.0015)|— 0.0012 |— 0.0036 0.0017 
wre S U i S U U U S lor 
Stability prediction | a S | U* = S| U,*=U » U, Us =S U; U; S an 
Actual Test poh aay buds «iy Sv Uy | ve | Uy Sv | Uv Uy NS ; 
i Fi A ae pla 
S = stable; U = unstable; U* = unstable unless (1) is large. for 
itn sp 
STABILIZING EFFECT OF COUPLED SURFACES ratio of times-to-damp-to-'/. amplitude cre 
eg . , : . short mode 
The purpose of a stabilizing surface on an airplane is ————— <} Th 
to restore equilibrium conditions following a disturb- long mode da 
ance, in a minimum of time and without perceptible Each of these modes will respond to a disturbance in ac- bet 
oscillations. The motion following a disturbance of a cordance with Eq. (12), as shown in Fig. 6. a 
one-degree-of-freedom vibrating system represented by The ratio of relative amplitudes of the Jong and short 
the equation ay + by + c = Ois: modes following a static displacement was found by 
W/vo = e7~*cos wot = e~"/"cos (1/8)(t Jt) (12) numerical solutions to be approximately as the ratio of 
the frequencies, and independent of the forward veloc- 
where +t = 2a/b represents a characteristic time;* 
£ = b/2aw represents the damping ratio as discussed ? 
in “Dynamic Stability of Coupled Surfaces’; and 
wo = Vo. ‘a represents the undamped natural frequency 
of the coupled tail, when operating on its own moment 
of inertia, m,x* only, as analyzed previously. 
In “Dynamic Stability of Coupled Surfaces’; the 
coupled tail was found to have two modes of motion, 
represented by short and long period oscillations fol- 
lowing a disturbance (cf. Table 1). 
; , wo short mode_ 2 
ratio of natural frequencies ————————_ > - 
wo long mode ~ 1 Fic. 6. Motion following a disturbance. 
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ity. These characteristics are believed to hold true in 
general because of the relative magnitude of terms re- 
sulting from “Optimum Proportions for Dynamic 
Stability,”’ which insures that the frequency separation 
of the modes shall be great and that both modes shall 
have adequate damping. This causes the short mode to 
vanish quickly, leaving the long mode predominant. 
Conversely, since the long mode is predominant, the 
procedure was correctly stated to increase the damping 
and natural frequency of the coupled tail, and to insure 
that the long mode should predominate. 
The damping ratio he for the long mode is 


— CraQSE Cra(p/2) SE a 
(U?/k*)m x? x" i 


(see Table 1). The damping ratio is seen to be inde- 
pendent of the velocity. The natural frequency, how- 


ever, 
Crag St K 
m,x* H 

increases linearly with the velocity, since g = (p/ 2)U*. 

If the coupled tail is operated as described previously 
(i.e., without the remainder of the airplane), in an air 
stream of velocity U, and the rudder is subjected to a 
sinusoidal forcing function cf frequency wy, the response 
of each mode is as shown in Fig. 7, point A. When the 
same tail is installed on an airplane, the natural fre- 
quencies wo of the two modes remain the same but the 
natural frequency of the airplane becomes less than the 
long mode of the tail alone, because of increased inertia 
and other effects. The natural frequency of the air- 
plane-plus-tail becomes the new forcing frequency wy 
for the tail, and as this frequency decreases with re- 
spect to the lowest wo of the tail alone, the ratio wy/wo de- 
creases and the rudder action varies as shown in Fig. 7. 
The true stabilizing effect, measured by the time to 
damp to '/; amplitude, then corresponds to point B, 
between the value determined for the dynamic condi- 
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Fic. 7. Response of coupled tail to forced oscillation. 


tion as in “Dynamic Stability of Coupled Surfaces,’ 
and the static conditions, as in “Static Stability of 
Coupled Fin and Rudder,’ depending on the change in 
frequency from “tail alone’ to ‘‘airplane-plus-tail.”’ 
For extremely long period motion, the static stability 
takes charge and the improvement factor becomes 
{1 + 7r(8/y)], corresponding to point C. 

The above qualitative discussion shows the general 
trend of the dynamic stability. However, in order to 
determine accurately the dynamic behavior of airplane- 
plus-tail, the equations of motion of the airplane, in- 
cluding the effect of the coupled tail, should be solved. 
The general solutions for (1) longitudinal motion and 
(2) lateral and directional motion, provide subjects for 
further analysis and research. 

The foregoing investigations assume no friction or 
free play in the control system. Friction and free play 
in the system wil] cause the rudder to lag behind the 
position it would have if these effects were absent. The 
result of this lag is a residua] moment acting on the sur- 
face at all times and easily recognized whenever the 
lift passes through the zero lift condition. This lag is 
contradictory to the original hypothesis that 8 is pro- 
portional to y. Furthermore, the residual moment due 
to lag is an unstable effect, tending to extract energy 
from the air stream and create an undamped oscilla- 
tion of the control, the amplitude of which increases as 
the lag angle increases. The effects of friction and free 
play may be corrected by reducing the linkage ratio, 
at the expense of control effectiveness [1 + 1(8/y)]. 
The quantitative effect of these factors remains a sub- 
ject for further investigation. 


SUMMARY 


In order to attain satisfactory results with the 
coupled tail, the surface geometry and airfoil contours, 
mass balance terms, coupling linkage, etc., must be in 
the correct proportions. Normal methods of construc- 
tion are considered to be satisfactory in general, pro- 
viding the design requirements for the coupled tail are 
met. 

The adaptation of manual control to the surface may 
be done in several ways. The major requirement is to 
avoid excess friction and free play in the system. The 
same precautions are necessary in developing the coup- 
ling linkage mechanism, which should be housed within 
the fuselage if possible (for aerodynamic cleanness). 

From the weight standpoint, the following comparison 
may be made: 

(1) Weight increase due to additiona] hinges and 
coupling linkage. 

(2) Weight decrease due to smaller surface and 
lower span. 

Although no actual weight calculations have been 
made, it is believed that, for a given stabilizing effect, 
the coupled tail may be produced at a lower weight than 


a normal tail. 
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The procedure for obtaining maximum stabilizing 
effect is: 

(1) Design for maximum static stability as in 
“Static Stability of Coupled Fin and Rudder.”’ 

(2) Check the damping criteria as in ‘“Summary of 
Damping Criteria.” 

(3) Maximize the damping effect by the procedure 
of ‘Optimum Proportions for Dynamic Stability,” by 
mass balances of fin and rudder, and proper selection of 
surface contours, linkage ratio, etc. 

(4) Low friction and free play are necessary for 
satisfactory results. 

By use of the above methods, it is possible to design 
stabilizing surfaces that are many times more effective 
than fixed surfaces of equal area. For example, the fol- 
lowing results may be expected from a good design: 


coupled tail - (: + 78) > 8 
fixed tail ¥ 


(cf. “Static Stability of Coupled Fin and Rudder’’) 


ratio of static stability 


coupled tail _ 
fixed tail 
'/,(H/K) > 6 


ratio of time to damp to '/, amplitude 


(cf. Table 1) 


When compared to a normal free rudder, with its 
downstream floating tendency, these ratios will be 
somewhat larger. The coupling linkage, in fact, causes 
a considerable upstream floating tendency. 

The advantages of coupled surfaces are shown by 
these comparisons. Whereas a normal surface may give 
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static stability in the free condition, the coupled surface, 
with the additional margin of static stability and in. 
creased damping, will restore equilibrium conditions 
much more quickly after displacement. It will-behave 
essentially like an automatic stabilizer. For this reason, 
properly designed coupled surfaces may be called 
“automatic control surfaces.”’ 
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Airplane Maneuvering Loads 


E. C. POSNER* 
Lockheed Aircraft Corporation 


ABSTRACT 


The tend in Army, Navy, and C.A.A. structural loading re- 
quirements is leading more and more toward rationalization. 
This involves the solution of the equations of motion to determine 
the dynamic response to certain arbitrary control surface move- 
ments. For example, the Army maneuvering tail load criterion 
depends on the airplane’s response to a given control movement 
applied in a fraction of a second, held there for a certain time 
interval, and then reversed to some percentage of the maximum 
value. The time interval involved is determined so that the 
maximum load factor reached is equal to the airplane limit factor. 
Likewise, the new rolling criterion consists of the application in a 
given time interval of an aileron angle of such a magnitude that 
a specified helical angle of roll is attained during the resulting air- 
plane response. The proposed C.A.A. vertical surface loading 
requirement consists of a rapid neutralization of the rudder angle 
at some rationally determined angle of sideslip. 

All these maneuvers have certain aspects in common: namely, 
they consist of a sequence of linearly varying control movement, 
followed by a constant deflection, and require the solution of the 
airplane’s equations of motion to determine the resulting response. 
Generalized graphical methods of analysis are available, but the 
author has found that the analytical methods are in general easier 
toapply. The analytical solution of these equations is presented 
here in a manner that it is hoped will be understandable to those 
not well versed with the mathematical theory involved. Rigor 
may therefore be sacrificed for the sake of clarity. 


NOTATION 


The following notation is similar to that used in standardized 
works on dynamic stability, such as references 2 and 3. It is 
therefore considered outside of the scope of this paper to go into 
the development of these derivatives, except for the derivative 
M,, which has not been generally used and needs some clari- 
fication. 


= longitudinal body axis—positive forward 
lateral body axis—positive to right 
vertical body axis—positive down 
force along z-axis, lbs. 

moment about x-axis, lb. ft. 

{ = moment about y-axis, Ib. ft. 

= moment about z-axis, Ib. ft. 

= velocity along x-axis, ft. per sec. 
velocity along y-axis, ft. per sec. 
velocity along z-azis, ft. per sec. ° 
acceleration along z-axis, ft. per sec.” 
angle of bank, rad. 

angle of pitch, rad. 

angle of yaw, rad. 

angular velocity of roll, rad. per sec. 
angular velocity of pitch, rad. per sec. 
angular velocity of yaw, rad. per sec. 
= angle of sideslip = v/U, rad. 

angle of attack = w/U, rad. 

control surface deflection, rad. 

’ = weight of airplane, lbs. 
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= mass of airplane = W’/g, slugs 

= air density, slugs per cu.ft. 

wing area, sq.ft. 

wing span, ft. 

tail area, sq.ft. 

tail length, ft. 

time, sec. 

m/pSb = mass factor 

m/pSU = time factor, sec. 

nondimensional velocity of roll = gy 

nondimensional velocity of pitch = 6v 

nondimensional velocity of yaw = yy 

t/vy = nondimensional time 

radius of gyration, ft. 

= angle of downwash, rad. 

= tail efficiency = ratio of dynamic pressure at tail to 
free stream 

An = load factor increment 

AP; = maneuvering tail load increment 


Hr z 
oeununw 
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Subscripts 
e = elevator 
r = rudder 
a = aileron 
t = tail 
m = maximum 


A“ THE ANGLES, velocities, and forces of the above 
notations are incremental values from the bal- 
anced condition at the start of the time history. An 
exception to this is the forward velocity, U, which is 
taken equal to the air speed at the start of the maneuver. 


Zu, Ly, Ls, Mv, Mz My My,N, N; = 
dZ dL aL 


, a] ’ e 
dw’ d¢’ dé 


Sn = LZu/pUS = Zyv/m 
= Myb/pUSk,? = Myv?U/mk,2u 


My = 
m;, = M;/pSk,? = M,vU/mk,? 
m, = M,/pUSk,? = M,v/mk,,? 


ms, = M,b/pU?Sk,? = M,v?/mk,?pu 
= N;b/pUSk, = N,?U/mk,2u 
n, =N,/pUSk,? = N,v/mk? 

= N,b/pU?Sk,2 = Nyw?/mk?u 
l, = L,/pUSk,? = L,v/mk,? 
L,b/pU?Sk,? = Ly?/mk,*u 


2 = 
oa 
| | 


oa" 
II 


Derivation of M;: 


M, “6 —'/pU?S1,C,, é 
M,, = dM/dw = (dM/da,)(da,/dw) 


de l; ($:)(=)(2) 
Aa; = t, a> Be spe — 
dt U,\da/\dw/\ dt 


But 
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where ¢, is the time for the downwash from the wing to 
reach the tail. Then 


dea,/dw = (l,/U,)(de/da)(1/U) 
and 
dM/da, = —'/»U?Sd,(dC,,/da;) 
Hence, 


Mz = —'/eV nSd2(dCz,/da,) (de/de) 


MANEUVERING HoRIZONTAL TAIL LOADS 


In the following derivation of the horizontal maneu- 
vering tail loads, the development is accompanied by 
generalized formulas so that the solution can be made 
applicable to dynamic problems of a more involved 
nature. 

As described in the abstract, the maneuvering hori- 
zontal tail loads depend on the dynamic response of the 
airplane to certain elevator movements. It is therefore 
necessary to set up the pertinent differential equations of 
motion and integrate these before the tail loads can be 
calculated. 

The response of the airplane to elevator deflections 
occurs in the longitudinal plane of symmetry and in- 
volves three degrees of freedom, x, 2, and 6. The cor- 
responding velocities are U, w and 6. However, the 
time interval during the part of the maneuver of interest 
is short so that the forward speed, U, can be regarded 
as constant. There then remain only two degrees of 
freedom and, hence, two equations of motion. 

Motion along z-axis: 


force = mass X acceleration 
wl» = m(w — 6U) (1) 


Motion about y-axis: 
moment = moment of inertia X angular acceleration 
5,M, + wM, + wM;, + 6M, = mk,?(d6/dt) (2) 


The term 6U represents the centrifugal acceleration 
due to the angular velocity of pitch, 6, and acts along 
the z-axis. 

It is convenient to convert into nondimensional form 
by multiplying the first equation by v/mU and the 
second by v?/mk,?. The resulting equations be- 
come 


a2. = (da/dr) — q (3) 
pdm, + pam, + m;(da/dr) + gm, = dq/dr (4) 


These are now placed in operator form by substitut- 
ing d to represent the operation d/dr. The equations 
then become 


(d — %)a—qg=0 (5) 
—(md + pmy)a + (d — m,)q = umd, (6) 


where 6, is a function of 7. 





1946 


The solution of these two simultaneous differentia] 
equations is arrived at by operating on each by a factor 
that enables elimination of one variable. Operating 
on Eq. (5) by d — m, and Eq. (6) by 1, then adding, 


da + (—2_ — mM, — m;,)da + (Zy9m, — uMy)a = 
HM 56, (7) 
This procedure is identical to that which would be 
followed if d were an algebraic multiplier instead of an 
operator. It is therefore possible to generalize by ob- 
taining Eq. (7) directly by means of determinants, 
The general formula for solution of Eqs. (5) and (6) is: 


a = [fi(d)/F(d)]6.(ry (8 
and 
q = [ fo(d)/F(d)] 6.(r) (9) 


where F(d) is the determinant made up of the operators 
of the variables on tue left-hand side of Eqs. (5) 
and (6): 


d — Sy —1 


= q? _ ae au 
—m,d — um, d— m, WBF eri, 


m;,)d + (2M, — wmy) 








fi(d) is the determinant with the right-hand coefficients 
substituted for the operators of a: 

0 —1 
um, d—m, 








= pM, 


f2(d) is the determinant with the right-hand coefficients 
substituted for the operators of q: 


d — Sy 0 


= um;(d — 
—md Pe UM um, b 3( Zw) 








The expansion of the general determinant 








Yo Me © 
Bi B, B; 
ie & 6] 


is 
A,B,C; + BiC,A3 + C,B3A2 — 
A;B,C, — B3C2A; — C3BiA, 


The solution for a@ is, hence, 





d? + (— 2 — m, — m;,)d + (Zw IMg — My) 
and for q is 


um,(d sons 2») 5(r) (11) 
ad? + (—Zy — mM, — m;,)d +(Z~mM, — My) 





The interpretation of Eqs. (10) and (11) requires 
explanation. It is first necessary to split up into partial 





fractions. This is done by the following expansion: 
Let 
d A B 
fd) _ Oe A 
Fd) d-m d-m d-ds 





H 





ifferential 
y a factor 
)perating 
idding, 


ax 
504 (7) 


would be 
sad of an 
e by ob- 
minants, 
1d (6) is: 


(8 


(9) 
perators 
Eqs. (5) 


#2 HM) 


fficients 


fficients 


C3BiA: 


(10) 


(11) 
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where Ai, Ae, As, ... . are roots of F(d) = 0 so that 
F(d) = (d — ¥)(d — A2)(d — As)... 


Cross-multiplying, 


fd) = A(d — x) (d—s).... + B@ — (A —).. 


, £ ~ 1 = wS.. 











But 
F’(d) = (d — de)(d — dz)... + 
(d — \y)(d — dAs).... + (d — A)(d — re)... 
and 
F’(\y) = (d — de) (d — dg)... 
F’ (Az) = (d- A) (d Rss. 
F’(A3) = (d — ¥i)(d — Ao)... 
Also 
f(m) = A(d — d2)(d — Ag)... 
f(\2) = Bid — )(d — As). . 
f(\s) = Cid — Mi)(d — a)... 
Hence, 
_ ee oe ee 
F’(\)’ F(X.) F’(Xs)" 
Then 


f(@)/F(d) = ONSO)/F’A)d — 2) (12) 


This is known as Heaviside’s expansion theorem. It 
should be noted that it breaks down when two roots of 
F(d) = 0 are equal, in which case other methods of ap- 
proach must be used. 

The integration of each term of this expansion is ob- 
tained as follows for each phase of the elevator move- 
ment. The first phase is assumed to consist of a linear 
variation of elevator angle with time until the time, 70, 
is reached, while the second phase consists of a constant 
deflection starting from time, 7. 


Linear Elevator Deflection, 5(r) = Kr 


Let 
u = 1r/(d — x) 
(du/dr) — \¥u = T 
(d/dr)(e~"u) = e~™r 
eulg = foe dr = — (e~™/d2)(A7 + 1)]5 
e~™u = (1/x*) — (e~™/d*)(Ar + 1) 
u = (e* — Ar — 1)/d? 
Hence, 
na Kr = Ky oe — dr —1) (13) 
and 


eo" —dr—-1 oO — hr — 1 


| | (14) 
12241 + B) * 922 + B) 





a= Kym) 


where F(d) = d? + Bd + C, F’(d) = 2d + B, and 
fi) = ym. 

It should be noted that the dependent variable u is 
assumed zero when the independent variable r is zero. 
Where the boundary conditions of the problem are such 


that this is not true, then Eq. (13) does not apply. 
Constant aeeaper Deflection Starting from Tine, 
To, (7) = Kro 


This phase of the motion can be considered to be 
made up of the difference between the linear motions 
Krand K(r — 7). Hence, 


f(d)/F(d)[Kr — K(r — 1)] = 














J(A) e? 
—rAr-1j — 
Ky) XFA) | lf 
A(r—ro) __ Xr vo — oh f(A) 
" pl Vt. BPO) 
[e"(1 — e~*”) — Aro] (15) 
and 
Air — To 
= é (1 =—*<¢ )—- AiTo 
° wm Qn +B) 
Aer —-— 2 nany 
é (J é ) ure (16) 
z?(2A2 + B) 


When the roots A, to A, are real, the development of 
Eqs. (13) and (15) is straight forward. However if, a 
root is complex, there is another root that is its con- 
jugate and the solutions must be converted into trigo- 
nometric form. This is done by converting all the com- 
plex terms into exponential form, combining, and then 
converting into trigonometric form. This can best be 
illustrated by going back to the original problem of 
Eq. (10). Here 


F(d) = d?+ Bd +C 


where B = —(2,. + m, + m;) and C = ym, — umy. 
If the roots of F(d) = O are complex quantities, let 


M= ath; wy =a-bi 
But 
\ = — (B/2) =V(B/2)* — 
Hence,a = — B/2andb = VC — a ora? +BP=C, 
and 


F’(d) = 24+B 








F’(A) = 2(a = bt) — 2a = +21 
Jud) = fil) = ums 
Hence, for the elevator motion 6(r) = Kr, from Eq 
(14), 
(a+bi)r (a + bi)r 7 
= K 
: uml (a + bi)*2bi 
(a—bi)r __ — * _ 
e (a . a ‘] (17) 
—(a — bi)?2bi 
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But the complex number a + bi expressed in polar form cos y = (e'? + e~")/2 
is sin y = (e? — e 7) /2t 
Va? + b? (cos y + isin 7) or cos y + isin y = e'’, and cos y — isin y = e~", 
Hence, 
where 


m a + bi = VCe” 
= =e = sin! —_ . 
y = tan~! (b/a) = sin (b/ VC) a —- bi = VGe7" 
Also the trigonometric functions can be expressed in (a + bi)? = Ce” 
exponential form as follows: (a — bt)? = Ce~"7 

















and 
ef"e" — V Cre” se ete * — V/Cre~” ae 
=Kym , 
patty - » { Ce" 2bi bs — Ce=*""2bi | 
"7 et Or—27) = ge ter—8y) P Pastas =! et? eh oe en 
wo) - aS)" 
bC 24 bV/C 2i bC2i 
= Kum ge sin (br — 2y) + eg ™ a . sin 2 
1 BC ma ee 
= Kym E sin (br — 2y) + + =| (18 
wl C*C ) 
For the elevator motion 6(r) = Kr — K(r — 7), from Eq. (16), 
(a+bdi)r — go (ardi)ro, __ : ° (a—bi)r sowie iey i fe . 
3 abel (l-—e ) (a + bi)to | ¢ (l—e é ) eed (19) 
(a + bi)? 2h1 —(a —b1)? 2hi 
()\(—— aan ——) ([\—— —s —— ( To \(- a. 
= Kums| { - : Ape : = | pal 
Cb 24 Cb 21 bVC 24 


sin [b(r — ro) — 2y] + a" 


ar a(r — 70) 


-. é 
= Kum \e sin (br — 2y) — 








a(r—ro) 
= Kum— [sin (br — 2y)(e“" — cos bro) + cos (br — 2y) sin bro] + Kum, 
= Kum, [(Ee*’~™ /bC) sin (br — 2y + &) + (70/C)] (20) 
where 


£ = tan—! [sin bro/(e"” — cos bro)] =sin ~' (sin b 1/E) 


and 








E= JV (e°” — cos bro)? + sin? bro = V er — 2e°" cos br> + 1 


The load factor increment An can be derived from the angle of attack increment a by the following relation: 


a = onset) () = t2)(S2)(D = - 4 
da/\W. 2\ da m g vg 
Knowing the angle of attack increment a as a function of nondimensional time r from Eqs. (14) and (16) if roots 


are real or Eqs. (18) and (20) if roots are complex, it is possible to construct a load factor versus time curve and 


choose the time increment Ar so that the limit load factor is just reached. 
The curve 0-1-2-4-5 of Fig. 1 is obtained by adding the respective load factor curves corresponding to the eleva- 


tor deflection 0-1-2-3 and the reversal motion 2-4-5. 
Once the time interval Ar is determined from the above plot, the tail loads can be easily obtained from the fol- 


lowing relation: 





AP, = n(dC,/da) i(p/2) U*S a; (22) 


*sin 2y = 2sin y cos y = 2ba/C 





Pe a a a 





(20) 


ation: 


(21) 


if roots 
rve and 


e eleva- 


the fol- 


(22) 
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= —| At | (3) 











An LIMIT L @ (5) 
An 
(0) 
O° 2 4 6 8 
i 9 
Fic. 1. Elevator and load factor time history. 
where 


da, de 
of— 1% eee 
” 4 * ( 3 “ 


and 6 = g/v = (da/vdr) — (az,/v) from Eq. (3). 














Elevator and tail load time history. 


Fic. 2. 





él, de \ (da\ {1; 
+ vet (EE) 2) 


The rate of change of angle of attack with time, da/dt = da/vdr, can be obtained by differentiation of the ex- 


pression for a. From Eq. (14). 


da/dr = Kum,y>N (e" — 1)/A(2d + B) (24) 


From Eq. (16) 


da/dr = Kum,>-Me™(1 — e~*)/A(2\ + B) (25) 


From Eq. (18) 


dr 
From Eq. (20) 
sa a(r—ro) 
ni Sime td Ge ~ 8-0 + eee 
dr bc 


g" 





cs = Kym, Ee cos (br — 2y) + = sin (br — 2y) + A = Kum, Ee. sin (br — y) + =| (26) 


cos (br — 2y + §)] = 


Ee*-™ 
bVC 





Kum, sin (b> — y + &) (27) 





The plot of Fig. 2 shows the tail load time history due 
to the elevator motion of Fig. 1. The tail load due to 
the respective components of Eq. (23) are shown. I is 
due to the airplane angle of attack corrected for down- 
wash at the tail; II is due to the angular velocity of 
pitch, or damping load; III is due to the elevator de- 
flection; and IV is due to the time lag of the down-wash 
in reaching the tail. 


It is evident that critical down tail load occurs at 
the start when the elevator is first fully deflected and 
the airplane has not had time to build up a positive 
angle of attack at the tail. The critical up load occurs 
when the elevator is reversed, since here, in addition to 
the load due down elevator, the angle of attack on the 
tail due to load factor is at its maximum and there is 
an appreciable damping load in the same direction. 
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MANEUVERING VERTICAL TaIL LOADS 


The maneuvering vertical tail loads are defined as the 
vertical surface loads due to the yawing motion of the 
airplane about the z-axis. If the banked altitude of 
the airplane is assumed to be zero at all times, as might 
be obtained by the proper aileron control on the part of 
the pilot, then there are only two degrees of freedom— 
the angle of sideslip, 8, and the angular velocity of 
yaw, y. 

Neglecting the side-force terms because these are of 
minor importance, the differential equations of motion 
in the yaw plane becomes: 

Motion along y-axis 


0 = m[(dv/dt) + YU] (28) 
Motion about z-axis 
6,N; + uN, + vN, sa mk,*(dy//dt) (29) 


Converting into nondimensional form by multiply- 
ing the first equation by v/mU and the second by 
v?/mk,*, these become, in operator form, analogous to 
Eqs. (5) and (6) 

dB+r=0 (30) 
—(yn,)B + (d — n,)r = pn,6, (31) 


where 4, is a function of r. 

It will be noted that these equations are similar to 
Eqs. (5) and (6) in that they are of second order with 
respect to 8, and fi(d) is — un; instead of um;. Hence, 
the solution for 8 and d8/dr can be written down di- 
rectly by comparing with the solutions for a and da/dr 
in Eqs. (14), (16), (18), and (20) and (24), (25), (26), 
and (27). 


Linear Rudder Deflection, 5(r) = Kr 


(a) Real roots 


Ae Ar — Xr — 1 
= —K > ek tell 32 
. 9 Luan MX + B) - 
dp ra e* an 
a fi ntact. 33 
de siti Pe + B) (33) 


(b) Complex roots 


p= —Kun, E sin (br — 29) +54 =| (34) 





d ar : 
- = —Kum| sin (br — y) + 3] (35) 


Constant Rudder Deflection Starting from Time, 7», 
5(r) = Kro , 


(a) Real roots 


(1 — e”) — Aro 
va > | oe 
b= Kuso 7. Han + B) oe 


dp ee aad @ ae e7 >”) 
Kum) AQ + B) = 
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(b) Complex roots 








a(r— 70) 
B= —Kyn| = sin (br — 2y + &) + “| (38) 
d E a(r—T0) : 
= = —Kun; We sin (br — y+) (39) 
where 
y = tan—'b/a = sin“ b/VC 
and 


sin bro : sin bro 
€ = tan~? ———_____ = sin? —_— 
e’” — cos bro 





where E = V e*” — 2¢"” cos bro + 1. 
The maneuvering tail load can be expressed in a 
manner similar to Eq. (22) 





= dC, p 
AP, = n( ) 2 U*S 81 (40) 
where 
dB, wl, 
= — 5, 41 
B, is, + e+ nt (41) 


and ¥ = (r/v) — (d8/vdr) from Eq. (30). 

It is desirable to derive an expression for the maxi- 
mum angle of sideslip developed in the oscillation re- 
sulting from a given rudder deflection. This is done by 
equatin Eq. (39) to zero and substituting the resulting 
value for r back in Eq. (38). 


btm — y+§ = ne 
T™ = (nx + y — &)/b (42) 


where m is the lowest positive integer that makes Tp 
greater than 7. Hence, 


2 (7m — 70) 
—Kun; (= sin (nx — y) + 4 


Il 


Bm 


be 
Eetm— 0) *) 
From Eq. (31) the value for the equilibrium angle 
Bo = —n;Kt0/n, (44) 
Hence, 
Bn inf Be (rm — 70) 4 Ee* (rm — 10) 
np Ge aan bo hima 1 Oe 
Bo T0 CVC + Cc nV C 


A typical time history of 8/8 is shown in Fig. 3 for 
a large conventional airplane. The locus of Bm/fp is 
also plotted for varying rudder deflection times, 79. It 
is interesting to note that the maximum angle of side- 
slip exceeds the equilibrium value by over 50 per cent 
for instantaneous rudder deflections and drops off 
slowly with increasing values of 7. It therefore ap- 
pears necessary to consider the time history in strength 
calculations even in the case of slowly applied rudder. 





—-> — oF je © - -. 





ed ina 


(40) 
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e maxi- 
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lone by 
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kes T, 
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Fic. 3. Slideslip time history. 


The neglect of side-force terms in Eq. (28) and the 
assumption that the angle of bank is zero at all times 
are conservative approximations as far as Fig. 3 is 
concerned, since the side force due to bank is restoring 
as well as the side force due to sideslip, both tending to 
reduce the maximum angle of sideslip. 


MANEUVERING ROLLING LOADS 


The maneuvering rolling loads are defined as the loads 
due to rolling motion about the x-axis. The rolling 
motion is produced by the application of a given aileron 
angle in a specified time, 7», and held constant thereafter. 
If the yawing moments induced by the aileron deflection 
and velocity of roll are considered baianced out by the 
rudder at all times, then there is only one degree of 
freedom involved, the angular velocity of roll, ¢. 

The corresponding differential equation in roll is: 


Motion about x-axis 
5,L5 + Ll, = mk,*(d¢/dt) (46) 


Converting into nondimensional form by multiplying 
by v?/mk,”, this becomes, in operator form, analogous 
to Eqs. (6) and (31). 


(d — l)p = ulsbg (47) 


where 6, is a function of r. 

By reference to the formulas already developed, the 
solution for p can be written down directly. Here, 
there is only one root of the complimentary equation 
and it therefore must be real. 


Linear Aileron Deflection, 5(r) = Kr : 
b = Kul,(e’”” — 1,r — 1)/l,? (48) 
dp/dr = Kul,(e’” — 1)/l, (49) 


Constant Aileron Deflection Starting From Time 7», 
5(r) = Kro 











pb = Kul, —— L ee (50) 
lpr — ga ste 
SS Rekorthermparen (51) 
T l, 
CONCLUSION 


It has been demonstrated that the maneuvering loads 
due to motions about all three axes of the airplane can be 
readily determined analytically for a sequence of con- 
trol movement varying linearly with time up to a speci- 
fied value and then held constant. The general for- 
mulas may be applied to any other type of motion that 
has a solution for the operation of Eq. (12). However, 
it will be found that, in general, the linear sequence is 
most adaptable to the limitations such as set up in the 
maneuvering horizontal tail load problem. 
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AIREX MANUFACTURING COMPANY, INC. 
ALLIS-CHALMERS MANUFACTURING COMPANY 
ALLISON DIVISION, GENERAL MOTORS CORPORATION 
ALUMINUM COMPANY OF AMERICA 
AMERICAN AIRLINES S 
AMERICAN OVERSEAS » INC. 
AMERICAN BOSCH CORPORATION 
AMERICAN PHENOLIC CORPORATION 
ASSOCIATED AVIATION UNDERWRITERS 
ASSOCIATED FOUNDRIES & MANUFACTURERS, INC. 
ATLAS SUPPLY COMPANY 
BAKER STEEL & TUBE COMPANY 
BEECH AIRCRAFT CORPORATION 
BELL AIRCRAFT CORPORATION 
BENDIX AVIATION CORPORATION 
BENDIX PRODUCTS DIVISION 
BENDIX RADIO CORPORATION 
ECLIPSE-PIONEER DIVISION 
EXPORT DIVISION 
FRIEZ INSTRUMENT DIVISION 
PACIFIC DIVISION 
SCINTILLA MAGNETO COMPANY 
THE BG CORPORATION 
BOEING AIRCRAFT COMPANY 
BOEING SCHOOL OF AERONAUTICS 
BREEZE CORPORATIONS, INC. 
ESSEX TOOL & DIE COMPANY 
FEDERAL LABORATORIES, INC. 
CAL-AERO TECHNICAL INSTITUTE, A DIVISION OF 
GRAND CENTRAL AIRPORT COMPANY 
CARRIER CORPORATION 
CESSNA AIRCRAFT COMPANY 
CHANDLER-EVANS CORPORATION 
CHASE NATIONAL BANK OF THE CITY OF NEW YORK 
THE CLEVELAND GRAPHITE BRONZE COMPANY 
THE CLEVELAND PNEUMATIC TOOL COMPANY 
AUTOMOTIVE-AIRCRAFT DIVISION 
CLIFFORD MANUFACTURING COMPANY 
CONSOLIDATED VULTEE AIRCRAFT CORPORATION 
NASHVILLE DIVISION 
STINSON AIRCRAFT DIVISION 
CONTINENTAL MOTORS CORPORATION 


CULVER AIRCRAFT CORPORATION 
CURTISS-WRIGHT CORPORATION 


THOMAS A. EDISON, INCORPORATED, INSTRUMENT 
DIVISION 








EDO AIRCRAFT CORPORATION 
THE ELECTRIC AUTO-LITE COMPANY 
ELECTROL INCORPORATED 
ENGINEERING AND RESEARCH CORPORATION 
ETHYL CORPORATION 
THE FAFNIR BEARING COMPANY 
FAIRCHILD CAMERA & INSTRUMENT CORPORATION 
FAIRCHILD AERIAL SURVEYS, INC. 
FAIRCHILD ENGINE AND AIRPLANE CORPORATION 
DURAMOLD DIVISION 
FAIRCHILD AIRCRAFT DIVISION 
RANGER AIRCRAFT ENGINES DIVISION 
FEDERAL TELEPHONE AND RADIO CORPORATION 
FIRESTONE AIRCRAFT COMPANY 
FLEETWINGS DIVISION, KAISER CARGO, INC. 
FLETCHER AVIATION CORPORATION 
G & A AIRCRAFT, INC. 
GENERAL AIRCRAFT EQUIPMENT, INC. 
GENERAL AVIATION EQUIPMENT COMPANY, INC. 
GENERAL ELECTRIC COMPANY 
GENERAL INSTRUMENT CORPORATION 
GENERAL MOTORS CORPORATION 
AC SPARK PLUG DIVISION 
AEROPRODUCTS DIVISION 
BUICK MOTOR DIVISION 
CADILLAC MOTOR CAR DIVISION 
CHEVROLET MOTOR DIVISION 
DELCO PRODUCTS DIVISION 
DELCO-REMY DIVISION 
EASTERN AIRCRAFT DIVISION 
FISHER BODY DIVISION 
FRIGIDAIRE DIVISION 
HARRISON RADIATOR DIVISION 
RESEARCH LABORATORIES DIVISION 
ROCHESTER PRODUCTS DIVISION 
THE B, F. GOODRICH COMPANY 
THE GOODYEAR TIRE & RUBBER COMPANY 
THE GRAY MANUFACTURING COMPANY 
GRUMMAN AIRCRAFT ENGINEERING CORPORATION 
GUIBERSON DIESEL ENGINE COMPANY 
W. & L. E. GURLEY 
HANNA ENGINEERING WORKS 
HASKELITE MANUFACTURING CORPORATION 
HAWAIIAN AIRLINES LIMITED 
THE HILLIARD CORPORATION 
THE INTERNATIONAL NICKEL COMPANY 
JACK & HEINTZ, INC. 
JOHNS-MANVILLE SALES CORPORATION 
KELLETT AIRCRAFT CORPORATION 
KENYON INSTRUMENT COMPANY, INC. 
WALTER KIDDE & COMPANY, INC. 
KOLLSMAN INSTRUMENT DIVISION, SQUARE D COM- 
PANY 


LANGLEY CORPORATION 

LAVELLE AIRCRAFT CORPORATION 

LEAR INCORPORATED 

LIBERTY AIRCRAFT PRODUCTS CORPORATION 


LINK AVIATION DEVICES, IN 
THE Li METER 


LONGINES-WITTNAUER WATCH COMPANY, INC. 


MAGNAFLUX CORPORATION 
METAL ee COMPANY 


THE: MARQUETTE 

THE GLENN L. MARTIN COMPAN 
THE W. L. MAXSON CORPORATION 
WARREN McARTHUR CORPORATION 


‘ASCO : 
MICROMATIC HONE CORPORATION . 
MINNEAPOLIS-HONEYWELL REGULATOR COMPANY 
MOORE DROP FORGING COMPANY 
NATIONAL CITY BANK OF NEW YORK 


NATIONAL CREDIT OFFICE, INC. 
THE NEW YORK AIR BRAKE COMPANY 
NORMA-HOFFMANN BEARINGS CORPORATION 
NORTH AMERICAN AVIATION, INC. 
NORTH AMERICAN AVIATION, INC., OF TEXAS 
NORTHROP AIRCRAFT, INC. 
NORTHWEST AIRLINES, INC. 
OWENS-CORNING FIBERGLAS CORPORATION 
PAN AMERICAN WORLD AIRWAYS SYSTEM 
THE PARKER APPLIANCE COMPANY 
PENNSYLVANIA-CENTRAL AIRLINES CORPORATION 
PESCO PRODUCTS COMPANY DIVISION, BORG-WARNER 
CORPORATION 
PHILLIPS PETROLEUM COMPANY 
PIONEER PARACHUTE COMPANY, INC. 
THE PURE OIL COMPANY 
REPUBLIC AVIATION CORPORATION 
J. P. RIDDLE COMPANY 
A. V. ROE CANADA LIMITED 
JOHN A. ROEBLING’S SONS COMPANY 
ROHR AIRCRAFT CORPORATION 
ROOSEVELT FIELD, INC. 
THE RYAN AERONAUTICAL COMPANY 
SCIAKY BROTHER: 
SCOTT AVIATION CORPORATION 
SHELL OIL COMPANY, INC. 
SIMMONDS AEROCESSORIES, INC. 
SKYDYNE, INC. 
SOCONY-VACUUM OIL COMPANY 
SOLAR AIRCRAFT COMPANY 
SPERRY GYROSCOPE COMPANY, INC. 
SQUARE D COMPANY 
STANDARD OIL COMPANY OF CALIFORNIA 
STANDARD OIL COMPANY (INDIANA) 
STANDARD OIL COMPANY OF NEW JERSEY 
SWEDLOW AEROPLASTICS CORPORATION 
THE TEXAS COMPANY 
THOMPSON PRODUCTS, INC. 
TINNERMAN PRODUCTS, INC. 
TITEFLEX, INC. 
TRANSCONTINENTAL & WESTERN AIR, INC. 
TRIPLETT & BARTON, INC. 
UNION CARBIDE AND CARBON CORPORATION 
BAKELITE CORPORATION 
HAYNES STELLITE COMPANY 
LINDE AIR PRODUCTS COMPANY 
NATIONAL CARBON COMPANY 
UNITED AIRCRAFT CORPORATION 
CHANCE VOUGHT AIRCRAFT DIVISION 
HAMILTON STANDARD PROPELLERS DIVISION 
PRATT & WHITNEY AIRCRAFT DIVISION 
PRATT & WHITNEY AIRCRAFT CORPORATION OF 
MISSO! 


URI 

SIKORSKY AIRCRAFT DIVISION 
UNITED AIR LINES, INC. 
UNITED STATES AVIATION UNDERWRITERS, INC. 
UNITED STATES RUBBER COMPANY 

DOMINION RUBBER COMPANY, LTD. 
THE VARIETY AIRCRAFT CORPORATION 
VICKERS, INC. 
VIDAL CORPORATION 
THE WACO AIRCRAFT COMPANY 
WARNER AIRCRAFT CORPORATION 
THE WEATHERHEAD COMPANY 
WESTERN AIR LINES, INC. 
WESTINGHOUSE ELECTRIC CORPORATION 
WESTON ELECTRICAL INSTRUMENT CORPORATION 
WYMAN-GORDON COMPANY 
YOUNG RADIATOR COMPANY 
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